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FOREWORD 


Work  described  in  this  report  was  accomplishes  by  the  Flight 
Mechanics  Division  of  the  Air  .  ^.e  fiiym.  Dynamic-  Laboratory  and  the 

Digital  Computation  Division  of  the  Aeronautical  Systems  Division  under 
Project  1451,  "Flight  Path  Analysis,"  Task  143109,  "Trajectory  and 
Motion  Analysis  of  Flight  Vehicles."  The  formulation  and  interim 
documentation  were  completed  by  Major  Urban  H.  D.  Lynch.  Programming 
was  accomplished  by  Mr.  Fay  0.  Young  of  the  Digital  Computation  Division 
(ASVCP) ,  Computer  Science  Center,  Aeronautical  Systems  Division. 

This  report  was  prepared  by  Major  Lynch  and  Mr.  John  J.  Dueweke  of 
the  High  Speed  Aero  Performance  Branch  (FXG),  and  combines  the  applicable 
portions  of  FDL-TDR-64- 1  ,  Part  I,  Volume  I,  with  the  interim  documen ta t i o 
Tile  overall  report  is  divided  into  four  parts: 


Part  I. 

Part  1  I  . 
Part  I  I  I 
Part  IV. 


Capabilities  of  the  Takeoff  and  Landing  Analysis 
Computer  Program 

Problem  Formulation 

User's  Manual 

Programmer's  Manual 


This  report  was  submitted  by  the  authors  in  June  1972. 


Th i  technical  report  has  been  reviewed  and  is  approved. 


PHILIP  P.  ANTONATOS 

Chief,  Flight  Mechanic:..  Division 

Air  Force  Flight  Dyne'  ics  Laboratory 
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ABSTRACT 


TOLA  is  an  acronym  for  T_akeO_ff  and  Landing  Analysis  digital  computer 
program.  This  part  describes  the  use  of  the  program. 

The  basic  program  provides  six  rigid-body  degrees  of  freedom  of 
aerospace  vehicle  motion  over  a  flat  planet  and  determines  the  response 
of  the  aircraft  to  control  inputs.  The  dynamics  of  up  to  five  independent 
oleo-type  struts  are  included  for  simulation  of  symmetrical  and  non- 
symmetrical  landings  as  well  as  drop  tests. 

A  maneuver  logic  is  programmed  to  provide  vehicle  guidance  in  the 
various  phases  of  the  problem;  it  determines  the  desired  trim  and 
position  in  the  glide  slope  and  provides  synthesis  and  attempted  com¬ 
pletion  of  necessary  flare  dynamics  for  a  prescribed  touchdown  velocity 
vector.  The  landing  roll  includes  wheel  spinup  and  braking,  thrust 
reversing,  spoiler  deployment,  and  system  failure  options.  The  takeoff 
roll  consists  of  acceleration  to  takeoff  speed,  followed  by  rotation  to 
takeoff  angle  of  attack. 

The  autopilots  attempt  to  obtain  smooth  response  with  no  overshoot 
and  provide  values  for  pitch,  yaw,  and  roll  control  surface  deflections, 
throttle  settings  for  one,  two,  three,  or  four  engines,  and  controlled 
(or  failed,  locked,  or  constant)  braking.  The  control  response  logic 
simulates  linear  control  system  rates  and  their  corresponding  lags  to 
des i r ed  va I ues . 

The  TOLA  Drogram  is  ideal  for  dynamic  tradeoff  studies  in  aircraft 
design,  landing  gear  design,  and  landing  techniques.  The  formulation  is 
programmed  for  the  CDC  6000  and  Cyber  70  Computer  Systems.  The  program 
is  programmed  in  Fortran  Extended  using  the  Scope  3-^  operating  system. 
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SECTION  I 
INTRODUCTION 

The  purpose  of  this  report  is  to  summarize  and  complete  the  docu¬ 
mentation  of  Project  1^+31 09-002  “Take  Off  and  Landing  Analysis',  (TOLA). 
TOLA  is  a  FORTRAN  modification  to  Option  2  (SDF-9)  of  FDL-TDR-6A- 1 , 

“S i x-Degree-of-Freedom  Flight  Path  Study  Generalized  Computer  program'', 
and  allows  comprehensive,  quantitative  calculation  of  aircraft  takeoff 
and  landing  performance.  Specifically,  this  report  shows  how  to  use  the 
TOLA  computer  program. 

The  reader  should  at  least  familiarize  himself  with  the  formulation 
documentation  before  attempting  to  use  the  TOLA  simulation.  This 
documentation  i^  contained  in  Volume  I  cf  this  report,  and  in  AFFDL- 
TR— 68- 111,  whn.h  formulates  the  equations  of  motion  for  a  series  of 
nonrigid  bodies. 

The  performance  analyst  will  need  to  cbtain  data  on  toe  various 
components  of  a  given  vehicle.  The  types  of  data  required  are  summarized 
under  TOLA  Inputs.  Data  Format  details  how  the  data  is  used. 
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SECTION  II 

TOLA  INPUTS  (Specify  Units) 


1.  VEHICLE  DATA 

Either  TOGW  or  landing  weight 

Ixx,  lyy,  Izz,  1X2  eg  location  (W.L.,  body  station) 

2.  LANDING  GEAR  DATA 
For  Each  Gear: 


Moment  of  inertia  of  a_  tire, 
wheel  ,  and  anything  else 
constrained  to  rotoie  with  the 
tire  about  the  o«!e  ,  (  Ik  )  MOMENT 


Mass 

strut 

outer 


of  all  portions  of 
which  move  relatlva  to 
sleeve  of  strut,  (mk)  MASS 


1 


WS5>W*s»¥l 


AFfDL-TR-7l- 155 
Part  II! 


v 


i; 

)' 


f 

f 


If  there  is  a  metering  pin,  A^  *■  A^(s)  implying  C  -  C(s). 


■ 

Coefficient  of  eliding 
friction  Between  piston  ond 
sleeve  of  strut,  Mu* 


Tire  Force  -  Defloctlon 
Curve  (  s  ) 


£K- 


F  -  a, 

i 


sb| 


from 

log  F  -  log  a(  +  b(  lo9  8 
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3.  AERODYNAMIC  DATA 

Aerodynamic  Reference  Area  AREFF 

Longitudinal  Reference  Length  (usually  M.A.C.)  D1RFF 
Lateral  Reference  Length  (usually  span)  D2RFF 

"Wind"  Axes  (with  and  without  spoilers) 

Co  vs.  Ot  for  various  elevator  deflections 

C,  vs.  ^  for  various  elevator  deflections 
“Body"  axes 

Cm  vs,,  oi  for  various  elevator  deflections 

(Origin  of  moment  axis  system  [W.L.,  B.S.j) 
vs_.  6  "aileron"  deflection 

Cn  vs^  6  for  various  rudder  deflections 

C  vs .  8 

y  - 

Cy  C(1|d  *  cnr 

Of  the  above  listed  coefficients,  the  first  3  are  mandatory.  The  next  3 
are  required  if  crosswinds  or  engine  failures  are  to  be  studied.  The 
rate  coefficients  may  be  input,  also,  if  available. 

Provision  is  made  for  coefficients  with  no  ground  effect  and  including 
full  ground  ef  fee  . 

A  Cp  due  to  gear. 

A.  PROPULSION  DATA 

"Installed"  chrust  as  a  function  ot  Mach  No.  and  "throttle  setting." 
Thrust  vectorls)  assumea  parallel  to  body  x-axis;  thus,  the  body 
y-  and  z-axis  displacements  of  the  thrust  vector('),  YN,  ZN,  are 
required.  Mass  is  assumed  constant  (therefore.  No  fuel  is  required). 
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5.  RUNWAY  DATA 

Length  (Rj  -  RLT 

Altitude  above  sea  level  -  RWHGR 

Elevation  angle  relative  to  horizontal  (Er)  -  ERDEG 

6.  OPTIONAL  DATA 


Reverse  thrust 

Drag  chute;  requires  chute  C^,  SR£.p,  and  location  of  chute  attachment 
point  relative  to  c.g. 

Wheel  braking,  including  "controlled"  braking  based  on  a  commanded 
"skid" 


Engine  fai lure(s) 


Control  variable  rate  lags 

Effects  of  rotating  machiner,  including: 

Rotation  rate  of  machinery  about  its  own  shaft  (u>r)  -  0MGRT 

Pitch  angle  of  shaft  relative  to  body  x-axis  in 

plane  parallel  to  body  x-z  plane  (9r)  -  THTRD 

Moment  of  inertial  of  rotating  machinery  about 

its  own  shaft  (l  )  AIXR7S 


7.  OTHER  REQUIRED  DATA 


The  lower  and/or  upper  limits  on  the  following  variables  are  required. 


Tire  deflection  (for  blowout) 

Braking  moments 

"$  skid"  for  controlled  braking 
Control  angular  acceleration  of  wheel  to  the 
angular  speed  required  to  maintain  a 
given  skid" 

Angle  of  attack 
Angle  of  attack  rate 


Computer  Symbol 

DELTAM 
MBL ,  MBU 


5 


0MECD1 

ALPHDS  ,  ALPHDL 
ALPDL 
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Type  Data 

Comt uter  Symbol 

Elevator  deflection 

DELQL, 

DELQU 

Aileron  deflection 

DELPL , 

DELPU 

Rudder  deflection 

DELRL , 

DELRU 

Elevator  deflection  for  takeoff,  land 

ing  roll 

DELQTO 

Elevator  deflection  rate  for  landing 

rol  1 

DELFDI 

Thrust:  Throttle  setting  above  which 

reverse 

throttle  should  not  be  engaged 

NB 

Time  to  activate  spoilers  (tsp) 

TSP 

Time  to  reverse  engines  (trv) 

TRV 

Time  to  release  drag  chute  ( tC(1 ) 

TCH 

Time  to  initiate  wheel  braking  (t^) 

TBK 

Linear  rates  for  simulating  control  response  lags: 

Elevator  rate 

DELHS 

Ai leron  rate 

DELA 

Rudder  rate 

DELRRD 

Throttle  rate 

NED1 

6 
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SECTION  I  I  I 
DATA  FORMAT 


1.  CARO  FORMAT 

The  program  input  routine  (READ)  expects  the  following  format. 


Card  Columns 

1  -  6 

7 

8-10 

1 1 

12-66 

67  -  72 

Field 

1 

1  1 

1  i  1 

IV 

V 

VI 

73  “  80 
VI  I  I 


Card  Field  I  -  Contains  the  symbolic  name  of  the  variable  which  data 


contained  in  Field  V  begins  loading.  Example:  Card  Column 


1 

GAM7D 

SIG7D 


12 

-1  .23 

90. 


Card  Field  II  Not  used 

Card  Field  III  -  Contains  the  words  DEC,  0CT,  BCD,  TRA,  INT,  or  is 
blank,  depending  on  the  type  of  data  to  be  loaded.  The  word  0CT  i nd i cates 
that  the  data  is  to  be  interpreted  as  octal  numbers.  The  word  BCD 
specifies  that  N  binary  coded  decimal  words  (N  punched  in  column  12) 
beginning  in  column  13  are  to  be  loaded.  The  word  TRA  denotes  to  the 
input  routine  that  all  data  has  been  input  and  to  return  control  to  the 
calling  program.  The  word  DEC  and  blank  are  equivalent  and  specifies 
that  data  loaded  is  decimal  data. 


^CT  Example 

Card  Column 

1 

8 

12 

NSMAIN 

OCT 

17 

BCD  Example 
Card  Column 

1 

8 

12 

REM 

BCD 

3SDF2-GEAR-MOD 

The  3  in  Column  12  specifies  3  words  wnere  each  word  Is  considered  to  be 
6  characters  including  blanks.  The  largest  number  of  6  character  words 
that  can  be  loaded  from  one  card  is  9.  Analysts  should  be  very  careful 
to  see  that  the  BCD  information  does  not  get  punched  Into  Field  VI.  This 
will  cause  an  input  error. 


'  -tw  ■  I  rnrr  ' 
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DEC  Example 

Card  Column  !  8 

VT  ABO  I  DEC  2,0.  ,  1 .67,20000. , I .67 

Note  that  the  first  character  in  Column  12  is  an  integer  and  the  input 
routine  will  load  only  one  integer  per  DEC  card  and  that  has  to  be  the 
first  number  punched  in  Field  V. 

VTAB01  DEC  2. ,0.  ,1 .67,20000. , 1 .67 

If  the  above  card  is  punched,  the  two  will  now  be  loaded  into  the 
machine  as  a  binary  floating  point  number.  Likewise,  the  other  numbers 
will  be  loaded  the  same  with  the  decimal  point  assumed  right  justified. 

If  anything  other  than  OCT,  BCD,  INT,  TRA  or  blank  appears  in  Field  II 
then  the  word  DEC  is  assumed. 


INT  Example 

Card  Column 

1 

8 

IP 

INT 

IP 

IP 

INT 

When  the  word  INT  is  "sed,  it  is  assumed  that  all  numbers  on  the  card 
will  be  loaded  as  integers.  If  only  one  integer  is  punched  per  card  the 
INT  may  be  punched  or  omitted. 

Card  Field  IV  -  Not  Used 

Card  Field  V 


The  actual  Input  data  tc  the  program  is  punched  In  the  Field  V.  DEC, 
INT  and  0CT  numbers  must  always  be  left  adjusted;  that  is,  data  must 
always  start  in  column  12  on  the  input  card.  All  numbers  are  separated 
by  a  "comma"  and  the  field  terminates  with  the  first  blank.  BCD  infor¬ 
mation  begins  in  Column  13  and  the  maximum  number  of  6-character  words 
per  card  is  nine.  Note  that  since  Field  V  ends  with  the  first  blank,  the 
user  may  punch  any  comments  in  the  remainder  of  the  field. 
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Card  Field  VI 


This  field  specifies  the  initial  subscript  of  the  data  in  Field  V. 

If  this  field  is  blank  an  initial  subscript  of  I  is  implied.  The  subscript 
may  appear  anywhere  within  the  field. 


Example 

Card  Column 

1 

12 

67 

PZEiyT 

30-470.4,41538.24,1)1538.21* 

1  or  blank) 

PZER0' 

42538.24,42538.24 

4 

In  the  example  above,  the  number  30470.4  is  loaded  into  the  first  cell  of 
the  array  PZERjGT.  On  the  second  card  42538.24  is  loaded  into  the  fourth 
cell  of  the  array.  The  one  and  four  punched  in  Field  VI  indicate  the 
subscript  for  the  array  PZER0. 

Card  Field  VI  I 


Not  used  as  far  as  the  input  routine  is  concerned.  This  may  be  used 
as  a  sequence  number  for  the  card. 

2.  TABLE  FORMAT 

The  various  types  of  tables  used  by  the  program  may  be  classed  as 
fol lows : 


One-D i mens i ona 1 

Tables 

Example  1 : 

NT  IRES 

nj  *  f (i  )  , 

Card  Column 

1 

12 

NSTRUT 

5 

NT  IRES 

4  .  ,6. ,6. ,6.  ,6 

1  ,2 . NSTRUT 


NSTRUT  ■  Fixed  point  number  which  is  the  number  of  struts  on  the 
aircraft.  For  example,  the  number  of  tires  on  strut  H2  is  6.  I.e., 
n2  “  f(2)  "  6. 

i  -  independent  variable  valies 
nj  Corresponding  dependent  variable  values 
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Example  2: 
Card  Column 


1 

INDAO! 
ATABO 1 


Aerodynamic  Data 

12 

1 

.0065,-00748 


INDA01  1  designates  that  there  are  data  in  ATAB01 .  The  first  data  point 
is  for  full  ground  effect;  the  second  data  point  is  for  no  ground  effect 
in  all  aerodynamic  tables. 


Two-Dimensional  Table 

Example:  VTAB0I 

Card  Column  1 

VTAB01 

N  =  Fixed  point  number  equal  to  2  times  the  number  of  independent 
variables.  For  a  20  point  table,  N  would  equal  20.  The  total  number  of 
machine  cells  required  for  this  table  is  41. 


XCG  =  f(M) 

12 

N »M]  ,XCG] ,M2  ,XCG2 ’  MN,XCG 


M|  =  Independent  variable  values 

xCg  “  Corresponding  dependent  variable  values 


N  -  Dimensional  Table 


Example: 
Card  Column 


T  =  F(N,Mn) 


I 

mow 

IT  1  OX 

TTAB10 

TTAB10 

TTAB10 


12 

NN 

NMN 

N| ,  N2  ,  Nj ,  . . . , 

MN, *  mn2.  mn3>  • 
^N] *  % (MM! 


nNN 

• • .  Mnnmn 

>  •  •  •  »  >^N  j 


TTAB10 


TN2’MN2’  TN2»mN2 . TnNN’Mn2 


67 

1 


A 


TTABI0 


Tn1’MNNMN’  Tn2,MNNMN’  TnNN’mNNMN 


NN  and  NMN  are  fixed  point  numbers  of  independent  variables.  Tw  ,Mli 

"1  N  ] 

Tf>  ,  M  are  values  of  independent  variables.  The  table  subscripts 
NN  nNNMN 


■^Displacement  numbers  required  when  any  table  exceeds  (1)  card.  See 
example,  page  9- 
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would  apply  to  the  N-d imens iona  I  table  as  well  as  the  two  dimensional 
The  total  number  of  machine  cells  required  for  an  N-d imens iona 1  table 
equals  NN  *  NMN  +  NN  +  NMN .  (Asterisk  here  indicates  multiplication.) 

Examples : 

C  r.  F (X,  Y)  NX  =  2  =  points  for  X 

NY  =  2  =  points  for  Y 

Machine  cells  required  2x2+2+2=8cel!s 

C  =  F (X,  Y,Z) 

Machine  cells  required 


NX  =  20  =  points  for  X 

NY  =  10  =  points  for  Y 

N2  =  1 5  =*  points  for  Z 

20  X  10  X  15  +  20  +  10  +  15  =  30^5  cells 


AFFDL-TR-71-155 
Part  III 


SECTION  IV 
DATA  PREPARATION 


1.  DATA  PREPARATION  -  GENERAL 

Before  preparing  data  which  may  actually  affect  the  trajectory,  a 
certain  amount  of  data  must  first  be  considered,  including: 

a.  Table  size  data 

b.  Identification  data 

c.  Required  data  input  for  each  case 

d.  Integration  data 


These  are  pointed  out  in  the  sections  to  follow, 
a.  Table  Size  Data 

Table  sizes  are  read  into  the  program  by  data  cards  that  follow  an 
STCASE  TAB  card.  All  tables  are  assumed  dimensional  (l)  unless  otherwise 
specified  by  data  input. 


Example:  Define  tables  TTA810  (30), 
ATAB01  (2) ,  and  VTAB01  (5) . 


Column 


1 

STCASE 

TTAB01 

ATAB01 

VTAB01 


8 

TAB 


TRA 


12 

30 

2 

5 


(define  tables  and  Base  case) 


(END  tab  le  def in i t ion) 

Stage  1  data  follows  this 
IRA  card. 


Tables  can  only  be  defined  in  a  Base  Case.  The  sum  of  all  table  sizes 
may  not  exceed  600  without  a  program  modification  to  COMMON  block  COMMON/ 
TABDIR/  in  routines  AUXR2 ,  TLU,  HI  HO,  TFFS,  and  AERO. 
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i 

j 

j 


i 


! 


* 


It 

1 


$ 


i 

! 


The  table  size  definition  in  the  Base  Case  must  be  sufficient  for 
any  table  in  a  succeeding  Merge  Case  with  a  larger  table  size  than  that 
found  In  the  Base  Case. 

b .  identification  data 

a.  Remarks.  For  identification  purposes,  three  lines  of  60 
characters  each  may  be  caused  to  be  printed  at  each  major  stage. 

EXAMPLE:  THIS  IS  A  SAMPLE  TRAJECTORY 

IDENTIFICATION  IS  MADE  BY  REM 


Field  1 

Field  1  1  1 

Field  V 

Field  VI 

REM 

BCD 

3THISbISbAb SAMPLE 

REM 

BCD 

2bTRAJECT0RY 

4 

REM 

BCD 

2 1  DENT  1 F 1  CAT  1 

6 

REM 

BCD 

3<J>NbISbMAi)EbBYbREM 

8 

Here,  the  leading  number  in  Field  V  Is  the  number  of  6  character 
BCD  words  on  this  card  (including  blanks).* 

b.  Case  Number.  For  specific  identification  of  each  case,  the 
case  number  is  printed  at  the  top  of  each  page. 

EXAMPLE:  CASE  3-02A4 

Field  I  Field  I  I  I  Field  V 

NCASE  BCD  1 3 .02A4 

where,  the  leading  I  in  Field  III  denotes  that  one  6-character  field 
is  to  be  read.  Only  one  BCD  word  is  allowable  in  the  program  for  case 
number  identification. 


*This  example  could  be  put  on  one  card,  but  has  been  segmented  here 
because  of  page  limitations. 
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c.  Required  Data  Input  for  Each  Case 


QTY . 

UNITS 

SYMBOL 

PT 

NOM 

VALUES 

REMARKS 

STCASE 

NCASE 

BCD 

0 

0 

Initial  card  (see 
case  number 

also  Merge  Data) 

REM 

BCD 

B 1  anks 

Remarks  (maximum 

30  words) 

NSTAGE 

NO 

1 

Stage  number  for 

initiation  of  problem 

t 

SEC 

TIME 

YES 

0. 

Initial  trajectory  time 

r* 

if> 

X 

SEC 

TIMSX 

YES 

0. 

Trajectory  time  at  start  of  stage: 
i  ,e.  ,  t  »  t-t 

STAGE  SX 

tIT^X 

SEC 

TMAX 

YES 

0. 

Maximum  estimated 

flight  time 

At 

SEC 

DELTS 

YES 

.  1 

Autopilot  time  interval 

1 NDATM 

NO 

1 

0 

Atmosphere  title 
deleted . 

and  computations 

1 

Computes  p,  ,  T 

J 

A-  V  V 

PRINT 

YES 

Print  interval 

PRTMIN 

YES 

Minimum  print  interval 

Integration  Data 

Symbol  Used 
by  READ  Routine 

Math 

Notation 

Symbol  used 
by  Integration 
Rout i ne 

Nominal 

Va  1  ue 

Remarks 

1 VARBH 

1 VARBH 

0 

Use  variable  step, 

“  1  ,  Use  F i xed  Step 

TIME 

t 

X 

0. 

Time  to  begin  integ. 

DELTS 

At 

DX 

.  1 

Time  Internal  to  ini 

AMINER 

At  . 

m  i  n 

DXMIN 

.001 

Mini  mum  At 

AMAXER 

At 

max 

DXMAX 

10000. 

Maxi  mum  At 

REtERI 

RELEKI 

RFLER1 

.00007 

Re i  .  error  tol  #1 

RELER2 

RELER2 

RELER2 

.000005 

Re  1  .  error  to] .  #2 

PRTMIN 

PRTMIN 

0. 

Print  Min Imum 

TIMER 

TIMER 

30 

Time  remaining  (sec) 
before  exiting  to 

Plot  routine 
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2.  DATA  PREPARATION  -  SUBPROGRAMS 

a.  Required  Data  -  Option  2  (2S0F) 


QTY. 

m 

UNITS 

SYMBOL 

PT. 

NOM. 

VALUES 

REMARKS 

SLUGS 

AMASS 

YES 

0. 

Mass  of  body 

Y 

DEG 

CAM7D 

YES 

0. 

Elevation  flight-path  angle 

a 

DEG 

SIG7D 

YES 

0. 

Horizontal  flight-path  angle 

h 

FT 

HGC7F 

YES 

0. 

Geodet i c  al ti tude 

*9° 

FT 

XGZ7F 

YES 

0. 

Origin  for  longitudinal 
di splacement 

Y 

go 

FT 

YGZ7F 

YES 

0. 

Origin  for  horizontal 
d  i  splacement 

V 

g 

FT/SEC 

VG77F 

YES 

0. 

Ground  referenced  velocity 

X 

a 

FT 

XG77F 

YES 

0. 

Initial  longitudinal  displacement 

Y 

q 

FT 

YG77F 

YES 

0. 

Initial  horizontal  displacement 

DEG 

PHIBD 

YES 

0. 

Body  Euler  angle,  roll 

DEG 

PSIBD 

YES 

0. 

Body  Euler  angle,  yaw 

0 

DEG 

THTBD 

YES 

0. 

Body  Euler  angle,  pitch 

P 

RAD/SEC 

PI77R 

YES 

0, 

Inertial  roll  rate,  body  axis 

q 

RAD/SEC 

OI77R 

YES 

0, 

Inertial  yaw  rate,  body  axis 

r 

RAD/SEC 

RI77R 

YES 

0. 

Inertial  pitch  rate,  body  axis 

AFX 

IBS 

DLFXP 

YES 

0. 

Generalized  axial  force 

AFy 

LBS 

DLFYP 

YES 

0. 

Generalized  horizontal  force 

AFz 

LBS 

DLFZP 

YES 

0. 

Generalized  vertical  force 

Ai.T 

FTLBS 

DLLTF 

YES 

0. 

General i zed  rol 1 i ng 
moment,  body 

AMt 

FTLBS 

DLMTF 

YES 

0. 

Generalized  pitching  moment, 
body 

ANt 

FTLBS 

OLNTF 

YES 

0. 

Generalized  yawing  moment,  body 

INDAPC 

NO 

0 

0 

Delete  a,B  print 

1 

Print  a,B 

INDADD 

NO 

0 

1 

Compute  and  print  u,p 

INDFPA 

NO 

0 

0 

Delete  print  o<  > ,0 

1 

Print  y ,  o 

INDFPR 

NO 

0 

0 

Delete  y ,o  calculations 

I  Compute  and  print  y  ,o 
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Required  Data  -  Option  2  (2SDF)  (Contd) 


UNITS 

SYMBOL 

PT, 

NOM. 

VALUES 

REMARKS 

INDGCR 

NO 

0 

0 

1 

Delete  great  circle  range 

Compute  and  print  R^ 

INDPLA 

NO 

0 

0 

1 

* 

Delete  platform  computation 

Compute  platform  angles  r  th 
platform  torqued  to  loce' 
geocentr  i  c 

INDACM 

NO 

0 

0 

1 

2 

Delete  accelerometer  calculations 

Compute  and  print  Ay  Ay  A^ 

P  P  P 

platform  accelerometer 

Compute  and  print  ax,ay>a2 

body  components  of  inertial 
acce ie rat  ion 

LBS 

FXB7P 

YES 

0. 

Initial  values  of 
summation  of  forces  in 

LBS 

FYB/P 

YES 

0. 

body-axes  system  including 

LBS 

FZB7P 

YES 

0. 

body  component  of  weight 

INDGRT 

NO 

0 

0 

1 

Delete  gimbal  rotation 

Cal cu 1  a t i ons 

Compute  0p,  Lp.ip  -'pitch,  yaw, 

INDRMC 

NO 

0 

0 

1 

Delete  rotating  machinery 
compute  t i ons 

Include  -otating  machinery 
compute  1 1 ons 

R.P.M. 

OMGRT 

YES 

0  . 

Rotational  rate  of  machinery 
about  i  is  own  shaft 

DEG 

THTRD 

YES 

0. 

Initial  pitch  ar.^le  of  shaft 
perpendicular  to  body  x-y  plane 

RAD/SEC 

BTABOI 

T 1 

- 

Table  of  pitch  ra,  of  shaft 
as  function  of  stage  time 

SLUG-FT2 

AIXR7S 

YlS 

0. 

Moment  of  inertia  of  rotating 
machinery  about  its  shaft 

1 NDWGT 

NO 

0 

0 

Delete  print  of  wtiuht 

1  Prim  wl 
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Required  Data  -  Option  2  (2SDF)  (Contd) 


QTY. 

UNITS 

SYMBOL 

PT. 

NOM. 

VALUES 

REMARKS 

INOUIN 

NO 

0 

0 

Winds  are  set  to  zero 

1 

Wind  velocities  specified  as  f(t) 

2 

Wind  velocities  specified  as  f(h) 

#0 

3* 

Wind  velocities  specified  as  f(Rg) 

*gw 

FT/SEC 

WTA30I 

T 1 

- 

Table  of  wind  velocity,  plus 
from  south  to  north,  f(lNDWIN) 

V 

gw 

FT/SEC 

WTAB02 

T 1 

- 

Table  of  wind  velocity,  plus 
from  west  to  east,  f(INDWIN) 

i 

gw 

FT/SEC 

WTAB03 

T 1 

- 

Table  of  wind  velocity,  plus 

from  up  to  downward,  f(lNDWIN) 


NOTE:  ^requires  INDGCR  =  I 
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b.  Vehicle  Physical  Characteristics  Subprogram  (VPCS) 


OTY . 

UNITS 

SYMBOL 

PT. 

NOM. 

VALUES 

REMARKS 

Vehicle  physical  characteristics 
deleted 

INDVPC* 

NO 

0 

0 

0 

1 

Vehicle  physical  characteri st  i  cs 
subprogram  used 

X 

c.g. 

FT 

XCGBF 

YES 

0 

- 

Longitudinal  body  C.G.  position 

Ax 

c  •  9  • 

1 

XX 

FT 

SLUG- FT2 

DXCGF 

AIXX8S 

YES 

YES 

0 

0 

- 

x  -  X 

c.g.  c.g.R£F 

Moment  of  inertia  about  body 
x  ax  i  s 

1 

yy 

SLUG-FT2 

AIYYBS 

YES 

0 

- 

Moment  of  inertia  about  body 
y  axis 

i 

ZZ 

SLUG-FT2 

AI22BS 

YES 

0 

- 

Moment  of  inertia  about  body 
z  Axis 

1 

*y 

SLUG-FT2 

AIXYBS 

YES 

0 

- 

Product  of  Inertia,  body  axes 

i 

xz 

SLUG-FT2 

AIX2BS 

YES 

0 

- 

Product  of  inertia,  body  axes 

1 

yz 

SLUG-FT2 

AIY2BS 

YES 

0 

- 

Product  of  inertia,  body  axes 

XX 

slug-ft2/$ec 

AIXXSI 

YES 

0 

- 

Rate  of  change  of  moment  of 
inertia 

i 

yy 

slug-ft2/sec 

AIYYSI 

YES 

0 

- 

Rate  of  change  of  moment  of 
inertia 

i 

zz 

slug-ft2/sec 

AI22SI 

YES 

0 

- 

Rate  of  change  of  moment  of 
inertia 

i 

xy 

slug-ft2/sec 

AIXYSI 

YES 

0 

- 

Rate  of  change  of  moment  of 
inertia 

1 

xz 

SLUG-FTZ/SEC 

AIX2SI 

YES 

0 

- 

Rate  of  change  of  moment  of 
inertia 

i 

yz 

SLUG-FT2/SEC 

AIY2SI 

YES 

0 

- 

Rate  of  change  of  moment  of 
inertia 

JC 

y 

FT 

ALYJOF 

YES 

0 

- 

Characteristic  distance  for  jet 
damp i ng  force 

2 

FT 

AL2J0F 

YES 

0 

- 

Characteristic  distance  for  jet 
damping  force 

i  r 

FT 

ALL JDF 

YES 

0 

- 

Characteristic  distance  for  jet 
damping  tnrce 

i 

m 

FT 

ALMJDF 

YES 

0 

- 

Characteristic  distance  for  jet 
damping  moment 

n 

FT 

ALNJDF 

YES 

0 

- 

Characteristic  distance  for  jet 
damping  moment 
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QTY . 

UNITS 

SYMBOL 

PT. 

NOM. 

VALUES 

REMARKS 

( INOVPC) 

(1) 

Vehicle  physical  characteristics 
subprogram  used  to  compute 
vehicle  characteristics  as 
f (MASS ,  TIME,  or  C.G.)  .  If 

INDTFF  «  0,  then  INDTSO  must  be 
input  as  1  to  read  values  below 
(to  INDIDT).  Otherwise  use 

INOVPC  =  0,  input  data. 

1 

*C'9-REF 

FT 

XCGRF 

YES 

0 

- 

Reference  longitudinal  body 
c.g.  position 

X 

c.g. 

FT 

VTABOI 

T 1 

- 

Table  of  c.g.  as  f(MASS) 

1 

XX 

SLUG-FT2 

VTAB02 

T 1 

- 

Table  of  moment  of  inertia  about 
x  axis,  f(m) 

1 

yy 

SLUG-FT2 

VTAB03 

T 1 

- 

Table  of  moment  of  inertia  about 
y  axi s ,  f (n) 

1 

zz 

SLUG-FT2 

VTAB09 

T 1 

- 

Table  of  moment  of  inertia  about 
z  axis,  f(m) 

INDXZS 

NO 

0 

0 

XZ  is  not  a  plane  of  symmetry 

0 

1 

XZ  is  a  plane  of  symmetry 

1 

xy 

blUG-FT2 

VTAB05 

T 1 

- 

- 

Table  of  product  of  inertia  f(m) 

1 

y* 

SLUG-FT2 

VTAB07 

T 1 

- 

- 

Table  of  product  of  inertial  fim) 

( INDXZS) 

1 

XZ  is  a  plane  of  -ymmetry  1  - 

1-0  *y 
yr 

INOXYS 

NO 

0 

0 

XY  is  not  a  plane  of  symmetry 

1 

XY  is  a  plane  of  symmetry 

0 

0» 

SLUG-FT2 

VTAB06 

T  1 

- 

Table  of  product  of  inertia  f(m) 

( INDXYSl 

1 

XY  is  a  plane  of  symmetry  1  -  0 

INOJOP 

NO 

0 

0 

Jet  damping;  4  »  t  »  4  « 

n  y 
z 

1 

V  !m'  ?n’  V^z  are  conputed 
f rom  tables 

1 

y 

FT 

VTAB08 

T  1 

- 

Table  of  characteristic  distance 
for  damping  force,  ) 

K 

i 

FT 

VTAB09 

T  1 

- 

Table  of  characteristic  distance 
for  jet  damping  force,  f(x£  ) 
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QTY. 

UNITS 

SYMBOL 

PT. 

NOM, 

VALUES 

REMARKS 

FT 

VTAB10 

T1 

- 

Table  of  characteristic  distance 
for  jet  damping  moment,  f(x£  ) 

*m 

FT 

VTABi ) 

T 1 

- 

Table  of  characteristic  distance 
for  jet  damping  moment,  f(xc  ) 

i 

n 

FT 

VTAB 1 2 

T 1 

- 

Table  of  characteristic  distance 
for  jet  damping  moment,  f(*c  ) 

INDIOT 

NO 

0 

0 

Inertia  derivatives: 

i  -  1  -  i  -  i  -  i 

xx  yy  zz  xy  xz 

i  -  0 

yz 

1 

Inertia  derivatives  are  computed 
from  tables 

1 

i 

XX 

slug-ft2/sec 

VTAB 1 3 

T 1 

- 

Table  of  inertia  derivative, 
f(ts) 

i 

yy 

slug-ft2/sec 

VTAB 1 4 

T 1 

- 

Table  of  inertia  derivative, 
f«s> 

i 

zz 

SLUG-FT2/SEC 

VTAB 15 

T 1 

- 

Table  of  inertia  derivative, 
f(S) 

i 

xy 

SLUG-r~2/:EC 

VTAB 1 6 

T 1 

- 

Table  of  inertia  derivative, 
f(ts) 

i 

xz 

slug-ft2/sec 

VTAB 1 7 

T 1 

- 

Table  of  inertia  derivative, 
f(t.) 

i 

yz 

slug-ft2/sec 

VTAB 18 

T 1 

- 

Table  of  inertia  derivative, 
f(ts) 

el8 

FT 

EPSI8 

YES 

0 

Incremental  error  in  X 

c.g. 

£  19 

SLUG-FT2 

EPS  19 

YES 

0 

Incremental  error  In  1 

XX 

e20 

SLUG-FT2 

EPS20 

YES 

0 

Incremental  error  In  1 

yy 

e2l 

SLUG-FT2 

EPS2I 

YES 

0 

Incremental  error  in  1 

zz 

c22 

SLUG-FT2 

EPS22 

YES 

0 

Incremental  error  In  1 

xy 

e23 

SLUG-FT2 

EPS23 

YES 

0 

Incremental  error  In  1 

xz 

e2<* 

SLUG-FT2 

EPS2*4 

YES 

0 

Incremental  error  In  1 

yz 
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c.  Aerodynamics  Subprogram  (3ACS) 


QTY. 

UNITS 

SYMBOL 

£!• 

NON. 

VALUE 

REMARKS 

NDAER 

NO 

0 

0 

Delete  aerodynamics  calculations 

0 

1 

Compute  aerodynamics  for 
controlled  aircraft;  moderate 
variations  (options  1-5) 

a 

LE3 

AA77P 

YES 

0 

Axial  force  (body  axis) 

y 

L8S 

YA77P 

YES 

0 

Side  force  (body  axis) 

nf 

LBS 

ANA77P 

YES 

0 

Normal  force  (body  axis) 

i 

FT-LBS 

ALA77F 

YES 

0 

Moment  about  body  x  axis 

m 

FT-LBS 

AMA77F 

YES 

0 

Moment  about  body  y  axis 

n 

FT-LBS 

ANA27F 

YES 

0 

Moment  about  body  z  axis 
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QTY. 

UNITS 

SYMBOL 

PT. 

NOM. 

VALUES 

(INDAER) 

1 

1 

S 

FT2 

AREFF 

YES 

0 

d, 

FT 

OIRFF 

YES 

0 

d2 

FT 

D2RFF 

YES 

0 

E1 

- 

EPS1 

YES 

1  . 

£2 

- 

EPS2 

YES 

0 

t 

£3 

- 

EPS3 

YES 

1 . 

e4 

- 

EPS4 

YES 

0 

e5 

- 

EPS5 

YES 

1. 

e6 

- 

tPS6 

YES 

0 

C7 

EPS7 

YES 

I . 

£8 

- 

EPS8 

YES 

0 

e9 

- 

EPS9 

YES 

1. 

e10 

- 

EPS  10 

YES 

0 

£ll 

- 

EPSI  1 

YES 

1 . 

£ 1 2 

- 

EPS  1 2 

YES 

0 

h 

FT 

AMAXH 

YES 

1.E6 

max 


REHARKS 

Compute  Aerodynamics  for 
Controlled  Aircraft,  Moderate 
Variations 

Reference  area 

Reference  length,  longitudinal 

Reference  length,  lateral 

Error  multiplier  for  C^ 

Incremental  error  in  CN 

Error  multiplier  for  C^ 

Incremental  error  in  C. 

A 

Error  multiplier  for  C^ 
Incremental  error  in  C^ 

Error  multiplier  for  C ^ 
Incremental  error  In  C ^ 

Error  multiplier  for  C 

m 

Incremental  error  In  C 

m 

Error  multiplier  for  Cn 
Incremental  error  In  C 

n 

Aerodynamic  Cut-Off  Altitude: 

h  >  h  Aero.  ■  0 

max 

h  «s  h 

max 

h>  295,275:  Constant  aero. 
h£  295.275:  Compute  aero. 


Constant  Aero  (Data  must  be  submitted  In  body  axes  system) 


C  A 

- 

CAMNU 

YES 

0 

Axial  force  coefficient 

CN 

- 

CNIMNU 

YES 

0 

Normal  force  coefficient 

CY 

- 

CYMNU 

YES 

0 

Side  force  coefficient 

CL 

- 

CLHNU 

YES 

0 

Rolling  moment  coefficient 

Cm 

- 

CMMNU 

YES 

0 

Pitching  moment  coefficient 

r 

V 

n 

- 

CN2MNU 

YES 

0 

Yawing  moment  coefficient 
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QTY.  UNITS  SYMBOL  PT.  NOM.  VALUES  _ REMARKS 


The  following  tables  are  two-dimensional  and  are  read  as  functions  of 
Mach  number,  M^.  The  table  indicators  are  nominally  zero  but  must  be 


i  nput 

as  1 

if  the  table 

is  to 

be  used: 

\ 

PER 

OEG 

ATAB01 

T2 

- 

SCA/3ct 

ca2 

or 

PER 

DEG2 

ATAB02 

T2 

- 

9CA/9ct2 

A 

C. 

ab 

PER 

DEG 

ATAB03 

T2 

- 

3C./98 

A 

CA  2 

B2 

PER 

DEG2 

ATABOA 

T2 

- 

9Ca/332 

CA 

6q 

PER 

DEG 

ATAB05 

T2 

- 

9CA/36q 

CA6q2 

PER 

DEG2 

ATAB06 

T2 

- 

oCa/962 

cAa3 

PER 

DEG2 

ATAB07 

T2 

32CA/9a96 

PER 

DEG2 

ATAB08 

T2 

- 

32C/\/3a96( 

CA66q 

PER 

DEG2 

ATAB09 

T2 

- 

3ZC./9B95 

A  i 

CNo 

• 

ATAB10 

T2 

- 

CN  AT  a  - 

CNa 

PER 

DEG 

ATAB1 1 

T2 

- 

3Cn/9o 

CNa2 

PER 

DEG2 

ATAB12 

T2 

- 

9CN/3a2 

Cn3 

PER 

DEG 

ATAB13 

T2 

- 

3Cn/3B 

V 

PER 

DEG2 

ATA014 

T2 

- 

3Cn/3B2 

C% 

PER 

DEG 

ATAB1 b 

T2 

- 

9CN/36q 

Sq2 

PER 

DEG2 

ATAB16 

T2 

3VMJ 

AFFDL-TR-71-155 
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QTY. 

UNITS 

SYMBOL 

PT. 

NOM.  VALUES 

REMARKS 

Cnc6 

PER  DEG2 

ATAB17 

T2 

- 

32CN/3ct3g 

°Na6q 

PER  DEG2 

ATAB18 

T2 

- 

32C../3a36 

N  q 

C"65q 

PER  DEG2 

ATAB 1 9 

T2 

- 

32C.,/303<5 

N  q 

Cu  • 

IMa 

PER  RAD 

ATAB20 

T2 

- 

3CN/3(ad,/2Va) 

CNax 

PER  RAD 
PER  FT 

ATAB2I 

T2 

- 

32CN/3(ad,/2Va)3x 

CNq 

PER  RAD 

ATAB22 

T2 

- 

3CN/3(qd1/2Va) 

Sx 

PER  RAD 
PER  FT 

ATAB23 

T2 

- 

32Cm/3 (qd . /2V  ) 3x 
n  i  a  i 

Cyo 

- 

ATAB24 

T2 

- 

Cy  AT  a  -  0  -  0° 

Cva 

PER  DEG 

ATAB25 

T? 

- 

3Cy/3ot 

Cyoc2 

PER  DEG2 

ATAB26 

T2 

- 

3Cy/3ct2 

Cy8 

PER  DEG 

. ATAB27 

T2 

- 

3C  /30 
y 

V 

PER  DEG2 

, ATAB28 

T2 

- 

3C  /302 
y 

PER  DEG 

ATAB29 

T2 

- 

3C  / 3  6 
y  r 

Cy^r 

PER  DEG2 

/.TAB30 

T2 

- 

3C  / 3  6  2 
y  r 

CVaB 

PER  DEG2 

ATAB31 

T2 

- 

32C  /3a36 

y 

Cy^r 

PER  DEG2 

ATAB32 

T2 

- 

32C  /3a36 
y  r 

CyB5r 

PER  DEG2 

ATAB33 

T2 

- 

32C  /303S 
y  r 

Cy  ‘ 

ye 

PER  RAD 

ATAB31* 

T2 

- 

3C  /3 (0d,/2V  ) 
y  i  a 

CY8x 

PER  RAO 
PER  FT 

ATAB35 

T2 

- 

32Cy/3(6d2/2Va)3x( 

2k 


7T'r«"-  h.-,  v 
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QTY. 

UNITS 

SYMBOL 

PT. 

V 

PER  RAD 

ATAB36 

T2 

CVrx 

PER  RAD 
PER  FT 

ATAB37 

T2 

C£ 

o 

- 

ATAB38 

T2 

C*a 

PER  DEG 

ATAB39 

T2 

CV 

PER  DEG2 

ATABAO 

12 

PER  DEG 

ATABAl 

T2 

<V 

PER  DEG2 

ATABA2 

T2 

CS 

PER  DEG 

.ATABA3 

T2 

PER  DEG2 

.atabaa 

T2 

C«-ag 

PER  DEG2 

ATABA5 

T2 

C£“6p 

PER  DEG2 

ATABA6 

T2 

c«-g<5p 

PER  DEG2 

ATABA7 

T2 

C*P 

PER  RAD 

. ATABA8 

T2 

C*r 

PER  RAD 

ATABA9 

T2 

^rx 

PER  RAD 
PER  FT 

ATAB50 

T2 

°mo 

- 

.ATAB51 

T2 

Cma 

PER  DEG 

.ATAB52 

T2 

Cma2 

PER  DEG2 

. ATAB53 

T2 

Cm 

m0 

PER  DEG 

ATAB5A 

T2 

S2 

PER  DEG2 

ATAB55 

T2 

NOM.  VALUES  _ REMARKS 

3Cy/8(rd2/2Va) 

32Cy/3(rd  /2V  )3x 

z  a  c.g. 

Cj>  AT  a  =  $  =  0° 

3C£/3a 

3C£/3a2 

3C£/9B 

3C£/3B2 

3C£/36p 

3C^/362 

32C£/3ct3B 

32C^/3a36p 

32C./3B36 

K  P 

3C£/3(pd2/2Va) 

3C2/3(rd2/2Va) 
32C£/3(rd2/2Va)3  xc 

C  AT  a  -  B  -  0° 
m 

3C  /3a 
m 

3C  /3a2 
m 

3C  /38 
m 

3C  /3$2 
m 
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QTY . 

UNITS  SYMBOL 

PT. 

NOM.  VALUES 

REMARKS 

Sq 

PER  DEG  .ATAB56 

T2 

- 

3Cm/96q 

Sq2 

PER  DEG2  . ATAB57 

T2 

- 

9C  /962 
m  q 

Cma6 

PER  DEG2  ATAB58 

T2 

- 

92Cm/9a96 

PER  DEG2  ATAB59 

T2 

- 

92C  /9a96 
m  q 

S6q 

PER  DEG2  ATAB60 

T2 

- 

92C  /9B96 
m  q 

Cma 

PER  RAD  .ATAB61 

T2 

- 

9C  /9(ad./2V  ) 
m  1  a 

S 

X 

PER  RAD  ATAB62 
PER  FT 

T2 

- 

92C  /9(ad,/2V  )9x  „ 
m  1  a  c.g 

H 

PER  RAD  . ATAB63 

T2 

- 

9C  /9(qd./2V  ) 
m  l  a 

Sx 

PER  RAD  ATAB6I) 
PER  FT 

T2 

- 

92C  /9  (qd  /2V  ) 9x 
m  1  a  c.g 

S 

ATAB65 

T2 

- 

C  AT  a  =  B  «  0° 
n 

a 

PER  DEG  ATAB66 

T2 

- 

9C  /9a 
n 

\2 

PER  DEG2  ATAB67 

T2 

- 

3C  /9a2 
n 

s 

PER  DEG  . ATAB68 

T2 

- 

ac  /as 

n 

S2 

PER  DEG2  . ATAB69 

T2 

- 

9C  /9B2 
n 

Sr 

PER  DEG  .ATAB70 

T2 

- 

3C  /96 
n  r 

Cf16r 

PER  DEG2  .ATAB7I 

T2 

- 

9C  /962 
n  r 

C"aB 

PER  DEG2  ATAB72 

T2 

- 

92C  /9a9B 

P 

Cna6r 

PER  DEG2  ATAB73 

T2 

- 

92C  /9a96 

n  r 

CnB6r 

PER  DEG2  ATAB7'» 

T2 

- 

92C  /3B9S 
rt  r 

Cne 

PER  RAD  ATAB75 

T2 

- 

9C  /9(8d,/2V  ) 
n  L  a 
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QTY. 

UNITS 

SYMBOL 

PT. 

NOM. 

VALUES 

REMARKS 

C"8x 

PER  RAD 
PER  FT 

ATAB76 

T2 

- 

32C  /3  ((3cL/2V  )3x  n 
n  2  a  c.g. 

Cflr 

PER  RAD 

. ATAB77 

T2 

- 

3C  /3(rd  /2V  ) 

n  4  3 

*'nr 

rx 

PER  RAD 
PER  FT 

ATAB78 

T2 

- 

32C  /3 (rd,/2V  ) 3x  n 
n  2  a  c.g. 

INDA80 

NO 

0 

0 

Table  not  used 

0° 

1 

Table  used 

- 

.ATAB80 

T3 

- 

CA  AT  a  -  6  -  0° 
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d.  Thrust  Subprogram  (TFFS) 


QTY . 

UNITS 

SYMBOL 

PT. 

NOM. 

VALUES 

REMARKS 

INOTFF 

NO 

0 

0 

Thrust  subprogram  deleted. 

Required  for  flight  plan 

Programmer  1  in  which  case  no 
additional  thrust  subprogram 
data  should  be  input. 

3 

Ai  rbreathing  engine 

3 

INOTSO 

MO 

0 

0 

1 

Converts  engine  thrust  to 
necessary  axes  systems 

(INDTSO) 

0 

0 

Deletes  thrust  conversion  to 
body  axes 

0 

L, 

FT-LBS 

ALT77F 

YES 

0 

Engine  roll  moment,  body  axis 

*T 

FT-LBS 

AMT77F 

YES 

0 

Engine  pitch  moment,  body  axis 

nt 

FT-LBS 

ANT77F 

YES 

0 

Engine  yaw  moment,  body  axis 

Tx 

LBS 

TXB7P 

YES 

0 

Longitudinal  thrust  component, 
body  axis 

T 

y 

LBS 

TYB7P 

YES 

0 

Lateral  thrust  component,  body 
axi  s 

T 

_ z 

LBS 

TZB7P 

YES 

0 

Vertical  thrust  component, 
body  axis 

~  0 

ZN( In) 

FT 

ZN 

YES 

0 

>n(in) 

FT 

YN 

YES 

0 

-3 

(INOTFF) 

3 

Airbreathing  engine 

T 

LBS 

TTABIO 

T5 

- 

Engine  thrust  as  f ( ,  N) 

IT  I0W 

NO 

0 

5 

Number  of  values  of  N  of  TTABIO 

ITIOX 

NO 

0 

Number  of  values  of  M^  for  TTABIO 

N  ( 1 N) 

- 

N 

YES 

0 

Throttle  setting 
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3.  DATA  ^REPARATION  -  STAGING 
a.  Staging 

For  a  trajectory  computation  which  requires  a  change  in  program 
computation  and/or  a  change  in  data  for  computations,  staging  steps  are 
used . 


Stage  Test.  Termination  of  a  stage  is  accomplished  with  the  following 
input  data: 


Field  I 

AiNCRS 

STEST 

and/or 

DECRES 

STESTD 


Field  I  1 1  Field  V 

BCD  nBBBBB  CCCC 

X.X,  YYY. 


BCD 


TRA 


nBBBB 

X.X 


whe  re 


n  ■  number  of  BCD  words  which  are  symbolic  in  the  directory 

(maximum  of  4  parameters  which  cannot  be  exceeded  (AINCRS) 
and  4  which  cannot  be  less  than  (DECRES). 


BBBBB 

CCC 

X.X 

YYY. 


)■ 

)• 


BCD  word  symbols  of  variables  computed  in 
program 


values  of  variables 
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A  stage  requires  two  sets  of  data,  each  followed  by  a  TRA  card. 

Primary  data 

FIELD  I  FIELD  III  FIELD  V 

INDSTR  0  nominal 

1  Reset  the  stage  time  to 

zero  and  update  the  stage 
number  by  one 


INDSP 


-N 


TRA 


where 

INDSP  =  indicator  for  subprogram  to  be  added  or  changed  from 
previous  stage. 

N  =  fixed  point  integer  value  of  indicator  desired  for  stage 

and  must  be  input  as  negative  for  proper  initialization. 


Secondary  data 

Any  data  required  by  addt d  or  changed  subprograms,  any  data  to  be 
changed  from  preceding  stage  data,  or  any  changes  in  auxiliary  com¬ 
putations,  are  added  followed  by  stage  test  data  and  a  TRA  card. 


Final  Stage.  If  a  particular  stage  is  to  be  the  last 
a  problem,  a  final  stage  indicator  is  required  prior 
card. 


or  final  stage  of 
to  the  final  TRA 


FIELD  I 


FIELD  V 

1 


INDSTF 
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b.  D  i  rectory-Symbol s  Available  for  Staging 

I  ri  order  to  stage  on  a  given  variable,  the  symbol  for  this  variable 
must  be  contained  in  the  program  directory.  Below  is  a  list  of  variables 
presently  available  for  staging.  They  are  listed  alphabetically  by 
symbol  name.  For  a  complete  listing  of  the  directory,  consult  the 
program  listings.  Values  listed  as  "not  used"  are  not  calculated  in  this 


version  of  the  program. 


SYMBOL 

QUANTITY 

UNITS 

AA77P 

a 

lbs 

AIXXBS 

'xx 

slug-f t2 

AIXXS1 

i 

XX 

s  1  ug-f t2/sec 

AIXYBS 

1 

xy 

slug-ft2 

AIXYS1 

i 

xy 

s  1  ug-f t2/sec 

AIXZBS 

1 

xz 

s 1 ug-f t2 

AIXZSI 

i 

xz 

s 1 ug-f t2/sec 

AIYYBS 

1 

yy 

s 1 ug-f  t 

AIYYS1 

i 

yy 

s 1 ug-f t2/sec 

AIYZBS 

1 

yz 

s 1 ug-f t2 

A!  YZS1 

i 

yz 

s 1 ug-f t2/sec 

AIZZBS 

1 

zz 

s 1 ug- f t2 

AIZZS1 

i 

22 

s 1 ug-f t2/ sec 

ALAMTD 

at 

deg 

ALA77F 

£ 

ft- 1 bs 

ALB77F 

L 

ft- lbs 

AL 1  FTP 

L 

1  bs 

ALPHD 

u 

deg 

ALPHRI 

ix 

rad/sec 

REMARKS 

Axial  force,  body  axis. 

Moment  of  inertia,  body  x  axis. 

Rate  of  change  of  moment  of  inertia. 
Product  of  inertia,  body  axes 
Rate  of  change  of  product  of  inertia. 
Product  of  inertia  body  axes. 

Rate  of  change  of  product  of  inertia. 
Moment  of  inertia,  body  y  axis. 

Rate  of  change  of  moment  of  inertia. 
Product  of  inertia,  body  axes. 

Rate  of  change  of  product  of  inertia. 
Moment  of  inertia,  body  y  axis. 

Rate  of  change  of  moment  of  inertia. 
Thrust  swivel  angle.  (Not  used) 

Moment  about  body  x  axis,  aerodynamic 
Moment  about  body  x  axis,  total. 

Lift  force  wind  axis. 

Angle  of  attack. 

Rate  of  change  of  angle  of  attack 


AFFDL-TR- 
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-71-155 

SYMBOL 

QUANT  1 TY 

UNITS 

REMARKS 

ALPTD 

a 

T 

deg 

Total  angle  of  attack.  (Not  used) 

ALT77F 

LI 

ft- 1  bs 

Engine  roll  moment,  body  axis. 

AMACH 

mn 

- 

Mach  number. 

AMASFS 

mf 

s  1  ugs 

Fuel  mass  consumed.  (Not  used) 

AMASF1 

mf 

s  1  ugs/sec 

Rate  of  change  of  fuel  mass.  (Not  used) 

AMASS 

m 

slugs 

Mass 

AMASS  1 

m 

s 1 ugs/sec 

Rate  of  change  of  mass. 

AMA77F 

m 

ft- lbs 

Moment  about  body  y  axis,  aerodynamic. 

AMB77F 

M 

ft- lbs 

Moment  about  body  y  axis,  total. 

AMT77F 

MI 

ft-lbs 

Engine  pitch  moment,  body  axis. 

ANAB7G 

n 

a 

g 1  s 

Axial  load  factor,  body  axis.  (Not  used) 

ANAZ7F 

n 

ft- lbs 

Moment  about  body  z  axis,  aerodynamic. 

ANA77P 

nf 

lbs 

Normal  force,  body  axes. 

ANB77F 

N 

ft-lbs 

Moment  about  body  z  axis,  total. 

ANGAMG 

n  Y 

g's 

Vertical  load  factor,  wind  axis. 

(Not  used) 

ANPSIG 

% 

g's 

Side  load  factor,  body  axis.  (Not  used) 

ANSIGG 

n 

a 

g's 

Side  load  factor,  wind  axis.  (Not  used) 

ANTHTG 

ne 

g's 

Vertical  load  factor,  body  axis. 

(Not  used) 

ANT77F 

nt 

ft-lbs 

Engine  yaw  moment,  body  axis. 

ANUA7F 

V 

f t^/sec 

Kinematic  viscosity  of  atmosphere. 

ANVW7G 

n 

V 

9's 

Axial  load  factor,  wind  axis. 

(Not  used) 

AXP7F 

p 

f  t/sec2 

Platform  acce lerat ion  ,  x  direction. 

AX77F 

a 

X 

ft/sec2 

Body  component  of  inertial  acce ler at  ion . 
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SYMBOL 

QUANTITY 

UNITS 

REMARKS 

AYP7F 

AYp 

ft/ sec^ 

Platform  acceleration,  z  direction. 

AY77F 

a 

y 

ft/sec2 

Body  component  of  inertial  acceleration. 

AZP7F 

A, 

ZP 

ft/sec2 

Platform  acceleration,  z  direction 

AZ77D 

a 

z 

ft/sec2 

Body  component  of  inertial  acceleration. 

BA77D 

8A 

deg 

Bank  angle.  (Not  used) 

BETAD 

6 

deg 

Angle  of  sidesl ip. 

BETARI 

6 

rad/sec 

Rate  of  change  of  sideslip  angle. 

B 1 770 

B 

deg 

Equatorial  angle  between  geocentric 
and  inertial  coordinate  systems. 

(Not  used) 

BP77D 

B 

P 

deg 

Equatorial  angle  between  inertial 
and  platform  coordinates.  (Not  used) 

CA 

CA 

- 

Axial  force  coefficient. 

CALPG 

Ca 

rad/g 

Gain  factor  for  angle  of  attack. 

(Not  used) 

CBETG 

CB 

rad/ g 

Gain  factor  for  angle  of  sideslip. 

(Not  used) 

CD 

CD 

- 

Drag  coefficient. 

CGAMG 

C 

Y 

- 

Gain  factor  for  flight  path  angle. 

(Not  used) 

CL 

C. 

U 

- 

Lift  coef f i c ient . 

CM 

Cm 

- 

Pitching  moment  coefficient. 

CN 

CN 

- 

Normal  force  coefficient. 

CP 

CP 

- 

Pressure  coefficient.  (Not  used) 

CRM 

C£ 

- 

Rolling  moment  coeificient. 

CY 

cY 

- 

Side  force  coefficient,  wind  axis. 

CYA 

CY 

- 

Side  force  coefficient,  body  axis. 

CYM 

C 

n 

- 

Yawing  moment  coefficient. 
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SYMBOL 

QUANTITY 

UNITS 

REMARKS 

DELPD 

6 

P 

deg 

Control  deflection  to  induce  a 
moment  about  the  x  axis. 

DELQD 

6 

q 

deg 

Control  deflection  to  induce  a 
moment  about  the  y  axis. 

DELRD 

6 

r 

deg 

Control  deflection  to  induce  a  moment 
about  the  z  axis. 

DLFXP 

AFx 

lb 

Generalized  axial  force. 

DLFYP 

AF 

y 

lbs 

Generalized  horizontal  force. 

DLFZP 

AF 

z 

lbs 

Generalized  vertical  force. 

DLLTF 

alt 

ft-  lb 

Generalized  rolling  moment,  body  axis. 

DLMTF 

AMt 

ft-lb 

Generalized  pitching  moment,  body  axis 

DLNTF 

AN_ 

ft-  lb 

Generalised  yawing  moment,  body  axis 

DRAGP 

D 

lbs 

Drag  force,  wind  axis 

DTC2R 

A0C2 

rad 

Pitch  attitude  command  correction  due 
to  temperature  and  temperature  rate. 
(Not  used) 

DXCGF 

AX 

eg 

ft 

X  -X 
eg  cgref 

DYNPP 

q* 

lbs/ft2 

Dynamic  pressure. 

OYNPP1 

q* 

lbs/f t^sec 

Rate  of  change  of  dynamic  pressure. 

FXB7P 

F 

X 

1  bs 

Summation  of  force  components,  body  x 
axis. 

FXE7P 

Fx 

xe 

lbs 

Summation  of  force  components,  planet 
Xg  axis.  (Not  used) 

FYB7P 

F 

y 

lbs 

Summation  of  force  components,  body  y 
axis. 

FYE7P 

\ 

lbs 

Summation  of  force  components,  planet 
Yfi  axis.  (Not  Used) 

FZB7P 

r. 

lbs 

Summation  of  force  co>nponents, 
body  z  axis. 
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SYMBOL 

QUANTITY 

UNITS 

REMARKS 

FZE7P 

Fz 

e 

lbs 

Summation  of  force  components,  planet 
Ze  axis.  (Not  used) 

GAMDD 

yd 

deg 

Geodetic  elevation  flight  path  angle. 
(Not  used) 

GAM7D 

Y 

deg 

Elevation  flight  path  angle. 

GXB7F 

9x 

ft/sec2 

Giavity  component,  body  x  axis. 

GXE7F 

9X 

e 

ft/sec2 

Gravity  component,  earth  Xe  earth  Xe 
axis.  (Not  used) 

GXG7F 

9X 

9 

ft/sec2 

Gravity  component,  geocentric  horizon 
coordinate. 

GXI7f 

9x, 

ft/sec2 

Gravity  component,  inertial  X  axis. 

GYB7F 

9y 

ft/sec2 

Gravity  component,  body  y  axis. 

GYE7F 

9Ye 

ft/sec2 

Gravity  component,  earth  Y  axis. 

(Not  used) 

GYI7F 

9Y, 

ft/sec2 

Gravity  component,  inertial  Y  axis. 

GZB7F 

9Z 

ft/sec2 

Gravity  component,  body  z  axis. 

GZE7F 

9ze 

ft/sec2 

Gravity  component,  earth  Ze  axis. 

(Not  used) 

GZG7F 

ft/sec2 

Gravity  component,  geocentric 
coordi nate. 

GZI7F 

9z. 

ft/sec2 

Gravity  component,  inertial  Z  axis. 

HGC7F 

h 

ft 

Geodetic  altitude. 

OMXGR 

UXg 

rad/sec 

Xg  component  of  planet  rotation  rate. 

(Not  used) 

OMZGR 

U)Z 

rad/sec 

Zg  component  of  planet  rotation  rate. 

9 

(Not  used) 

PA77P 

P 

lbs/ft2 

Atmospheric  pressure. 

PH  1  AO 

*A 

deg 

Aerodynamic  roll  angle.  (Not  used) 

PHIBD 

<P 

deg 

Body  Euler  angle,  rc’I. 
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SYMBOL 

QUANTITY 

UNITS 

REMARKS 

PHIGD 

*9 

deg 

Geodetic  latitude.  (Not  used) 

PHILD 

*L 

deg 

Geocentric  latitude.  (Not  used) 

PHI  PD 

*P 

deg 

Angle  between  body  axis  and  platform 
axis . 

PHiTD 

deg 

Rotation  angle  of  plane  of  thrust 
swivel.  (Not  used) 

PI77R 

P 

rad/sec 

Inertial  angular  rate  of  body  about 
x  axis . 

PI77R1 

P 

rad/sec^ 

Rate  of  change  of  p. 

PSIBD 

'  4* 

deg 

Body  Euler  angle,  yaw. 

PSIPD 

V 

deg 

Angle  between  body  axis  and  platform 
axis . 

QA770 

qA 

deg/sec 

Aeroelastic  body-bending  angular 
rate.  (Not  used) 

QI77R 

q 

rad/sec 

Inertial  angular  rate  of  body  about 
y  axis . 

QI77R1 

q 

rad/ser.^ 

Rate  of  change  of  q. 

RA77D 

rA 

deg/sec 

Aeroelastic  body-bending  angular 
rate.  (Not  used) 

RD77N 

rd 

n.mi . 

Total  distance  traveled  over  planet 
surface.  (Not  used) 

RG77N 

R 

9 

n.ml . 

Great  circle  range. 

RHOAS 

P 

slugs/f 

Atmospheric  density. 

RI77R 

r 

rad/sec 

Inertial  angular  rate  of  body  about 
z  axis . 

RI77RI 

r 

2 

rad/sec 

Rate  of  change  of  r. 

RPHLF 

V 

ft 

Local  planet  radius.  (Not  used) 

R777F 

R 

ft 

Distance  from  center  of  planet  to 
body . 

SIDEP 

Y 

lb 

Side  force,  wind  axis. 

36 


AFFDL-TR-71-155 
Part  ill 


SYMBOL 

QUANTITY 

UNITS 

REMARKS 

SIGDD 

°D 

deg 

Geodetic  horizontal  flight-path 
angle.  (Not  used) 

SIG7D 

a 

deg 

Horizontal  flight  path  angle. 

SIG7R1 

a 

rad/sec 

Rate  of  change  of  horizontal  flight 
path  angle. 

TA77R 

T 

°R 

Temperature  of  the  atmosphere. 

TE77P 

T 

lbs 

Engine  thrust. 

THL7D 

0L 

deg 

Longitude.  (Not  used) 

THTBD 

Q 

deg 

Body  Euler  angle,  pitch. 

THTPO 

9 

P 

deg 

Angle  between  body  and  platform  axis. 

THTRD 

9 

r 

deg 

Pitch  angle  between  rotating  machinery 
axis  and  platform  axis. 

TIME 

t 

sec 

Flight  time. 

TIMES 

*s 

sec 

Stage  time. 

TSTGR 

Te 

°R 

Equilibrium  stagnation  temperature. 
(Not  used) 

TSTGRI 

U 

°R/sec 

Equilibrium  stagnation  temperature 
rate.  (Not  used) 

TS77R 

Ts 

°R 

Skin  temperature.  (Not  used) 

TS77R1 

fs 

° R/sec 

Rate  of  change  of  skin  temperature. 

(Not  used) 

TVACP 

Tvac 

1  bs 

Vacuum  thrust. 

TXA7P 

Tv 

XA 

lbs 

Longitudinal  thrust  component,  wind 
axis . 

TXB7P 

TX 

lbs 

Longitudinal  thrust  component,  body 
axis  . 

TXE/P 

TXp 

€ 

lbs 

Longitudinal  thrust  component, 
planet  Xe  axis.  (Not  used) 

TYA7P 

Ty 

ya 

lbs 

Lateral  thrust  component,  wind  axis. 

37 


AFfDt-TR-71-155 
Part  III 


SYMBOL 

QUANTITY 

UNITS 

REMARKS 

TYB7P 

ty 

lbs 

Lateral  thrust  component,  body  axis. 

TYE7P 

Tve 

lbs 

Lateral  thrust  component,  planet 

Yfi  axis.  (Not  used) 

TZA7P 

TZ 

A 

lbs 

Vertical  thrust  component,  wind  axis. 

TZB7P 

TZ 

lbs 

Vertical  thrust  component,  body  axis. 

TZE7P 

TZ 

e 

lbs 

Vertical  thrust  component,  planet 

Ze  axis.  (Not  used) 

U777F 

u 

ft/sec 

Inertial  velocity  component,  body 
x  axis. 

U777F1 

u 

f t/sec^ 

Rate  of  change  of  u  component  of 
velocity. 

VA77F 

VA 

ft/sec 

Ai  rspeed. 

VD77F 

VD 

ft/sec 

Velocity  decrement  due  to  drag. 

(Not  used) 

VGRVF 

Vgrav 

ft/sec 

Velocity  decrement  due  to  gravity. 
(Not  used) 

VG77F 

V9 

ft/sec 

Ground  referenced  speed. 

VI77F 

V 

ft/sec 

Inertial  velocity. 

VP77F 

V 

P 

ft/sec 

Velocity  decrement  due  to  rocket 
nozzle  back  pressure.  (Not  used) 

VS77F 

Vs 

ft/sec 

Local  speed  of  sound. 

VTHEF 

VTHEO 

ft/sec 

Theoretical  velocity  increment  due  to 
Tvac.  (Not  used) 

V777F 

V 

ft/sec 

Inertial  velocity  component,  body 
y  axis. 

V777FI 

V 

f  t/sec^ 

Rate  of  change  of  v  component  of 
velocity. 

WTR7P 

wt 

lb 

Weight  of  vehicle. 

W777F 

w 

ft/sec 

Inertial  velocity  component,  body 
z  >xi  s . 
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SYMBOL 

QUANTITY 

UNITS 

REMARKS 

W777FI 

w 

f  t/secZ 

Rate  of  change  of  w  component  of 
veloci ty. 

XCGBF 

XCg 

ft 

Longitudinal  body  center-of-gravi ty 
pos i tion. 

XD77D 

xd 

deg 

Great  circle  downrange.  (Not  used) 

XD77N 

XD 

n .mi . 

Great  circle  downrange.  (Not  used) 

XE77F 

Xe 

ft 

Planet  referenced  coordinate. 

(Not  used) 

XGW7FI 

X 

gw 

ft/sec 

North-south  wind  velocity. 

XG77F 

X9 

ft 

Local  geocentric  displacement. 

XI77F 

X 

ft 

Inertial  coordinate. 

YA77P 

y 

lbs 

Side  force,  body  axis. 

YD77D 

yd 

deg 

Great  circle  crossrange.  (Not  used) 

YD77N 

yd 

n.mi 

Great  circle  crossrange.  (Not  used) 

YE77F 

Ye 

ft 

Planet  referenced  coordinate. 

(Not  used) 

YGW7FI 

Ygw 

ft/sec 

East-west  wind  velocity. 

YG77F 

V9 

ft 

Local  geocentric  displacement. 

YI77F 

Y 

ft 

Inertial  coordinate. 

ZE77F 

Ze 

ft 

Planet  referenced  coordinate. 

(Not  used) 

ZGW7FI 

zgw 

ft/sec 

Vertical  wind  velocity. 

ZG77F 

Z9 

ft 

Local  geocentric  displacement. 

ZI77F 

Z 

ft 

Inertial  coordinate. 
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4.  DATA  PREPARATION  -  DATA  MERGING 

To  eliminate  reproducing  a  large  number  of  cards  for  successive  cases 
with  similar  input  data,  a  data  merge  of  a  "base"  case  and  a  succeeding 
case  or  cases  can  be  made.  A  "base"  case  is  a  complete  set  of  input  data. 
A  "merge"  case  is  data  which  differs  from  the  "base"  case.  Control  of 
this  merging  facility  is  controlled  by  the  "STCASE"  or  start  case  card. 

The  STCASE  Card 


The  first  card  of  each  case  must  be  the  STCASE  card.  This  card  has 
the  same  fields  as  the  genera!  card  format  and  has  the  following 


meanings : 

Field  1 

Field  II 

Field  1 1 1 

Meaning 

STCASE 

Execute  this  case  as  a  base  case. 

STCASE 

NOXEQ 

This  case  is  a  base  case,  but  no 
computation  of  it  is  desired. 

STCASE 

MRG 

This  case  is  to  be  merged  with  the 
last  base  case  which  has  preceded 

it. 

STCASE 

MRG 

NOXEQ 

This  is  a  merge  case,  but  is  not 
to  be  executed.  This  set-up  should 
be  used,  since  physically  removing 
data  cards  will  accomplish  the 
same  function. 

In  order  to  illustrate  the  use  of  the  data  merge,  a  sample  case 
could  be  set  up  as: 


TTeTcTT 

mr 


BASE  CASE 


First 

Stage 

Next 

Stage 


Las  t 
Stage 


(data) 

(data) 


INDSTF 


IT 


TRA 


TRA 


TRA 


I  I  I 


Field  I 
St CASE 


INDSTF 


MERGED  CASE 


I  I 


I  I  I 


MRG 

Input  differs  from 
base  case 

TRA 

I nput  d i f fers  f rcm 
base  case 

TRA 


TRA 
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Note  that  there  is  a  one-to-one  correspondence  of  TRA  cards  between 
the  base  and  its  merged  case.  The  cards  of  the  merged  case  would  follow 
directly  behind  the  base  case.  Succeeding  merge  cases  will  use  the  last 
preceding  base  case.  Any  combination  of  base  and  merge  cases,  regardless 
of  execution  or  no  execution,  may  be  run  on  a  given  job.  The  final  stage 
of  both  the  base  case  and  each  succeeding  merge  case  must  contain  the 
INOSTF  card.  A  merge  case  may  contain  any  number  of  stages  irrespective 
of  the  number  of  stages  in  the  base  case.  However,  there  must  be  a 
one-to-one  correspondence  of  stages  up  to  the  point  of  departure. 
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SECTION  V 
INPUT 

The  data  deck  setup  will  be  explained  in  the  order  that  the  deck 
actually  appears.  A  typical  data  deck  input  is  shown;  however,  depending 
on  the  kind  of  problem,  all  thvi  data  shown  here  may  not  be  required. 

This  will  be  explained  later. 

1.  BASIC  SOP- 2  DATA 


1 

8 

12 

$DATA 

STCASE 

TAB 

ATAB01 

2 

ATAB02 

2 

ATAB10 

2 

ATAB1 1 

2 

ATAB12 

2 

ATAB15 

2 

ATAB16 

2 

ATAB27 

2 

ATAB28 

2 

ATA8l*3 

2 

ATAB44 

2 

ATAB48 

2 

ATAB51 

2 

ATAB52 

2 

ATAB53 

2 

ATAB56 

2 

ATAB57 

2 

ATAB61 

2 

ATAB63 

2 

ATAB68 

2 

ATAB69 

2 

ATAB70 

2 

ATAB71 

2 

ATAB77 

2 

ATAB80 

2 

TTAB10 

29 

VTAB01 

5 

VTAB02 

5 

VTA BO 3 

5 

VTABOl* 

c 

✓ 

VTAB06 

5 

FTAB02 

5 

FTAB03 

50 

CTAB01 

17 

CTAB02 

17 

CTAB03 

17 

CTAB0i< 

17 

WTAB01 

1 1 

WTAB02 

1 1 

WTAB03 

1 1 

FTABOI 

TRA 

90 

AFFDL-TR-71-155 
Part  I  I  I 


The  STCASE  TAB  card  (note  columns)  indicates  the  beginning  of  the 
data  deck.  This  card  indicates  that  all  tables  used  by  the  program 
will  now  be  sized.  This  is  followed  by  cards  indicating  the  tables 
names  and  sizes.  There  are  basically  five  kinds  of  tables  that  must 
be  sized:  aerodynamics,  thrust,  vehicle  physical  properties,  wind 
tables,  and  landing  mod  tables.  There  are  80  possible  aerodynamics 
tables  named  ATAB01  -  ATAB80.  The  relationship  between  the  table 
name  and  a  particular  aerodynamic  coefficient  is  shown  on  pages  26-30. 
Only  those  tables  actually  used  need  to  be  sized.  Because  of  the 
assumptions  made  in  Reference  1  (pgs  1 46 ,  224-225),  each  aerodynamic 
table  has  a  size  of  two.  The  thrust  table  is  TTAB10.  Because  of  the 
assumptions  made  in  Reference  1,  thrust  table  has  only  two  independent 
variables  (i.e.,  throttle  setting  and  Mach  number).  There  are  seven 
possible  vehicle  physical  properties  tables  named  VTAB01  -  VTAB07  (see 
pgs  21-22).  There  are  three  landing  mod  tables  named  FTAB01,  FTAB02, 
FTAB03.  FTAB01  is  the  force-deflection  curve  (i.e.,  fR(6R),  P9  99, 
Reference  1)  for  a  single  tire.  There  is  an  option  to  use  an  equation 
for  the  force-deflection  curve  in  which  case  FTAB01  is  not  needed. 
FTAB02  is  the  runway  perturbation  table  (i.e.,  e(XRR) ,  pg  110, 

Reference  l).  FTAB03  is  the  tire-runway  coefficient  of  friction  table 
(i.e.M^,  pgs  99-I0C,  Reference  l).  The  CTABO 1 -CTAB04  tables  are  for  the 
strut  orifice  coefficients  (see  Paragraph  2,  Landing  Gear  Modification 
Data).  The  last  card  indicating  a  table  name  and  size  is  followed  by  a 
TRA  card  (see  pg  13).  The  table-size  cards  between  the  STCASE  TAB  and 
TRA  cards  can  have  any  order. 
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Column:  1  8  12 


STCASE 

RFM 

BCD 

3SDF2-GEAR-MOD 

NCASE 

BCD 

IG-MO!) 

IVARBH 

0 

TMAX 

100. 

AMINER 

.00005 

PRTM 1 N 

.02 

AMAXER 

.1 

DFLTS 

.  1 

PRINT 

.  1 

AMASS 

19762.852 

HGC7F 

120.824 

TH1  6D 

3.5896 

GAM7D 

-3. 

VG77F 

220. 

RWHGR 

0. 

INDAPC 

1 

INDADO 

1 

INDEPA 

1 

INDPLA 

1 

INDACM 

2 

INDGRT 

1 

INDWGT 

1 

INDVPC 

1 

XCGRF 

0. 

VTAB01 

2,0. ,1.67, 20000. ,1.67 

VTAB02 

2,0. ,31 . 8E+06, 20000. ,31 .8E+06 

VTAB03 

2 , 0 . , 29 • 4E*  06 , 20000 . , 29 . 4E+06 

VTAB04 

2,0. ,57- 0E+06, 20000. ,57-0E+06 

INDXZS 

1 

VTAB06 

2,0. ,2. 39F+06 , 20000 . , 2 . 39E+06 

,  REM,  and 

NCASE 

cards  are  discussed  on  pgs  14  and  15. 

:ven  cards 

from 

IVARBH  to  PRINT  contain  integration 

'ARBH  must 

have 

a  zero  value.  TMAX  is  the  maximum  time 

span  (in  seconds)  expected  for  a  particular  calculation.  If  the  actual 
time  variable  in  the  program  exceeds  TMAX,  the  problem  will  stop.  AMINER 
is  the  minimum  integration  interval  allowed  by  the  variable  step 
Runge-Kutta  (VSRK)  integration  technique.  It  is  wise  to  set  AMINER  to  a 
very  small  value  to  allow  the  VSRK  the  flexibility  to  pick  as  small  an 
integration  interval  as  it  needs.  There  is  no  need  to  "second  guess" 
a  reasonable  value  for  AMINFR  and  run  the  risk  of  picking  a  value  too  large. 
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The  data  on  the  next  four  cards  PRTMIN,  AMAXER,  DELTS,  PRINT  must-repeat 
mus t-conform  to  the  following  relationship: 

PRTMIN  <  AMAXER  <  DELTS  <  PRINT. 

PRINT  is  the  output  print  interval.  DELTS  is  the  time  interval  at  which 
all  lutopi lot  calculations  are  updated  (that  is,  so  far  as  the  VSRK 
intec, ration  technique  is  concerned,  all  control  variables  remain  constant 
for  :.he  time  interval  DELTS).  AMAXER  is  the  maximum  integration  interval 
allowed  by  the  VSRK  integration  technique.  PRTMIN  is  the  lower  bound  on 
the  output  print  interval.  Data  output  on  a  PRTMIN  condition  will  only 
occur  when  the  VSRK  has  picked  an  integration  interval  smaller  than 
PRTMIN  and  remains  at  this  integration  interval  long  enough  for  the  time 
sran  to  exceed  PRTMIN.  Output  on  a  PRTMIN  condition  is  an  indication 
that  the  data  is  changing  rapidly.  The  PRTMIN  option  was  added  primarily 
to  output  the  landing  gear  data  at  a  sufficient  frequency  to  see  the 
detail  gear  response.  If  PRTMIN  is  given  a  small  value  (i.e.,  0.01  sec, 
et:.)  one  must  be  willing  to  accept  a  large  volume  of  data  output.  When 
PR'MIN  is  made  zero,  data  output  will  only  occur  at  the  PRINT  interval. 

The  cards  beginning  with  AMASS  down  to  and  including  INDWGT  are 
equired  data  for  the  original  SDF-2  Option  (see  pages  1 7~ 1 9 ) .  All  input 
data  on  those  pages  is  applicable  except  the  following: 

1)  INDADD  must  have  the  value  1 

2)  INDPLA  must  have  the  value  1 

3)  INDACM  must  have  the  value  2 

A)  INDGRT  must  hove  the  value  1 

The  cards  beginning  with  INDVPC  and  ending  with  VTAB06  are  the 
vehicle  physical  properties  (see  pgs  20-23)-  INDVPC  must  have  the 
value  1.  Data  format  for  the  VTAB  tables  is  shown  on  pg  8.  Since 
Reference  1,  pg  )A6,  makes  the  assumption  of  constant  mass,  only  two 
data  points  are  needed  in  each  of  the  VTAB  tables.  RWHGR  is  an  added 
input  and  represents  runway  altitude  (ft) 
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Column:  1  12 


INDAER 

1 

AREFF 

6200. 

DIRFF 

30.9 

D2RFF 

219. 

INDA01 

1 

ATABO! 

.0065,-00748 

INDA02 

1 

ATABO 2 

•437E-3..542E-3 

INDAIO 

1 

ATABIO 

1.34,1.12 

INCA) 1 

1 

ATAB1 1 

.1112, .1038 

INDAI2 

1 

ATAB12 

-.224E-2,-. I263E 

INOAI5 

1 

ATABI5 

-.00746, -.008 19 

INDA1 6 

1 

ATABI6 

.000)42, .000158 

INDA27 

1 

ATAB27 

-.0236, -.0236 

INDA28 

1 

ATAB28 

.267E-3, -267E-3 

INDA43 

1 

ATAB43 

.230E-3, -230E-3 

INDA44 

1 

ATAB44 

.212E-4, .212E-4 

INDA48 

1 

ATAB48 

-.45, -.45 

INDA51 

1 

ATAB51 

.0439, .2439 

INOA52 

1 

ATAB52 

-.0312, -.0241 

INDA53 

1 

ATAB53 

.  630E-3 , • 293E-3 

INDA56 

1 

ATAB56 

.0315,-0346 

INDA57 

1 

ATAB57 

-.6E-3.-.667E-3 

INDA61 

1 

ATAB6I 

-8. ,-8. 

INDA63 

1 

ATAB63 

-24. ,-24. 

INOA68 

1 

ATAB68 

. 182E-2, . 182E-2 

IN0A69 

1 

ATAB69 

. 101E-3, ■ 10IE-3 

INDA70 

1 

ATAB70 

. 1 58E-2 , . 158E-2 
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Column:  1  12 


INDA71 

1 

ATAB71 

-.346E-5.-.346E-5 

INDA77 

1 

ATAB77 

-.323, -.323 

INCA80 

1 

ATAB80 

.139,-155 

The  next  series  of  cards  beginning  with  INDAER  and  ending  with 
ATAB80  contain  aerodynamic  data.  INDAER  may  have  one  of  two  values. 

A  zero  value  deletes  the  aerodynamics  calculations  (note  that  no 
aerodynamic  data  is  needed  here).  The  only  place  this  option  would  have 
value  is  in  the  study  of  vehicle  impact  where  there  is  no  atmosphere. 

A  one  value  for  INDAER  causes  the  aerodynamic  calculations  to  be  made. 
Pages  24-30  show  the  possible  aerodynamic  input  data.  Figure  8  of 
Reference  i,  specifically  the  block  entitled  "Summation  of  Aerodynamic 
Forces  and  Moment  Coefficients,"  shows  the  form  of  the  aerodynamic 
coefficient  equations.  Due  to  the  requirements  of  the  autopilot  (see 
Reference  1,  pgs  224-225),  input  data  for  the  C^  and  CN  coefficients 
must  be  the  CD  and  CL  coefficients,  respectively,  in  the  pseudo  wind 
axes  system.  The  Cy,  C^,  Cm,  Cn  coefficients  must  be  body  axis  data. 

As  discussed  in  Reference  1  (pg  225),  only  two  data  points  are  needed 
for  each  aerodynamic  table  used,  the  first  data  point  is  for  full 
ground  effect  and  the  second  is  for  no  ground  effect. 

Column:  1  12  67 


INDTFF 

3 

INDTSO 

1 

mow 

5 

1 T 1  OX 

4 

TTABIO 

-2. ,-l . ,0. ,  1 .  ,2. 

TTABIO 

0. , . 1 , . 2 , .  3 

6 

TTABIO 

-13000. ,-800. ,0. ,2200. ,37400. 

10 

TTABIO 

-20500. ,-3400. ,-900. ,1300. ,34600. 

15 

TTABIO 

-30500. ,-7500. ,-1700. ,500. ,3)600. 

20 

TTABIO 

-4300. ,-14300. ,-2600. ,-400. ,28400. 

25 

The  next  ser!es  of  cards  beginning  with  INDTFF  and  ending  with  the 
last  TTA610  cird  pertain  to  engine  thrust  data.  The  indicator  INDTFF 
must  (see  Section  IV  para  3b  for  exception)  have  the  value  3-  INDTSO 
should  always  have  a  one  value  (see  pgs  21  and  30-  IT10W  indicates  the 
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number  of  throttle  settings  stored  in  the  thrust  table.  Its  value  must 
be  at  least  5  (that  is,  one  must  store  thrust  data  at  least  at  the  five 
throttle  settings  -2.,  -1.,  0.,  !.,  2.  -  see  Reference  1,  pg  225).  1 T 1  OX 
indicates  the  number  of  Mach  data  points  in  the  thrust  table.  The  first 
TTAB10  card  lists  the  throttle  setting  points  in  ascending  order, 
followed  by  the  Mach  points  in  ascending  order.  The  rest  of  the  TTAB10 
table  contains  the  thrust  data  for  one  engine  for  fixed  Mach  (start  vith 
lowest  Mach)  and  ascending  throttle  setting.  The  last  data  in  the  table 
is  for  the  highest  Mach  setting  and  ascending  throttle  setting.  The 
"displacement  numbers"  on  the  far  right  side  of  the  TTAB10  cards  is 
explained  on  pg  $. 


2.  LANDING  GEAR  MODIFICATION  DATA 


NSTRUT  is  the  number  of  landing  gears  on  the  aircraft.  It  is 
synonymous  with  the  variable  K  in  Reference  1  (pg  74).  The  program  is 
presently  limited  to  a  maximum  value  of  5  for  NSTRUT.  A  minimum  of  3 
struts  is  necessary  to  support  the  aircraft  on  the  runway.  MASS  Is  an 
array  of  the  masses  (i.e.,  mk,  see  Reference  I,  pgs  71  and  84)  in  slugs, 
of  the  main  strut  (i.e.,  strut,  wheel  frame,  tires,  etc.)  of  each  landing 
gear.  In  each  gear  array  there  must  be  NSTRUT  data  points.  The  order  of 
the  data  points  is  arbitrary  but  must  be  consistent  in  all  ge;r  arrays. 

If  the  data  for  the  array  exceeds  that  which  can  be  put  on  one  card, 
displacement  numbers  must  be  put  on  the  following  cards  for  that  array. 
The  arrays  RX ,  RY ,  and  RZ  are  the  body  coordinates  R^^,  Rfc  ,  R^, 
respectively,  in  feet,  of  the  landing  gear  position  vector  (R^k)0 
(Reference  I,  pg  74  and  Eq  114,  pg  82)  from  the  vehicle  nomlna’  mass 
center.  It  is  these  three  gear  position  arrays  that  determine  the 
appropriate  order  of  the  data  in  the  other  landing  gear  arrays.  RX  is 
positive  in  the  forward  \ose  direction.  RY  is  positive  oui  the  right  wing 
direction.  RZ  is  positive  down.  THETAD  is  the  array  of  angles  Gk 
(Reference  1,  pgs  74  and  80)  in  degrees.  EFIDEG  is  the  runway  elevation 
angle  Ek  (Reference  1,  pg  8l)  1'  degrees.  RGR  is  the  fixed  distance 
RgR  (Reference  I,  pg  81)  in  ft.  NTIRE5  is  an  arrav  of  the  number  of 
tlrer  nk  (Reference  I,  pg  99/  on  each  strut  axle.  RZERO  is  an  array  of 
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Col umn : 


1 

12 

67 

NSTP.UT 

5 

MASS 

A6. 62 1A9, 132.0942, 132.0942, 132. 0942 

MASS 

132. 79A2 

5 

RX 

62. 3^333,-1.  ‘a.  ~i.A2,-19. 75333, -19. 75333 

RY 

0., -12. 9375, 12. 9375. -12. 9375, 12. 9375 

RZ 

0. ,0. ,0. ,0. ,0. 

THETAD 

0. ,0. ,0. ,0. ,0. 

ERDEG 

0. 

RGR 

1917. 

NT  IRES 

A. ,6. ,6. ,6. ,6. 

RZERO 

2.01 1666,2.01 1666, 2. 011  666,2. Jl 

RZ'-T 

2.01 1666 

5 

OF  .M 

.83333,. 83333,. 83333,. 83333,. 83333 

RLT 

1 . E+05 

IFD 

1 

Al 

1 25^*72. 2, 1^1992.9, 141992.9, 1  A1992.9 

Al 

141992.9 

5 

Bl 

1.32, 1.293, ..293, 1.293, 1.293 

FTAB02 

2,0. , 0 . , 1 . E+G  7,0. 

FTAB03 

1A,  . ,0. , .01 , ,i, .029, .3, .035,. A 

FTAB03 

. 065,. 5,. 09,. 6,. 15.. 76,. 2,. 8,. 2A,. 76 

10 

“TAB03 

.325, .65. .5, .6, .7, .56,1  53,100. ,.53 

20 

MOMENT 

13.5,15.0,15.0,15.0,15.0 

MB 

o 

d 

d 

0 

0 

RF 

18.1 3066, 10.55416, 18.55416, 18.>.  ;6 

RF 

18.55416 

5 

V  Z 

PZERO 

304/0.4,41538. 24,4,538.24 

PZE'.O 

A1538.2A,A1538.2A 

A 

VZERO 

-  5579861 , 1 . 162662, 1 . 162662 

VZERO 

1 . 162662,1 . 1 62/6 n 

A 

A 

.A92222, .722222, .722222 

A 

.722227, .722222 

P20 

3758AO. , 3  4C1 5A . 4 , 3A6 1 5A . A 

P20 

3A6 1 5A .  A ,  3A6 1  £s4 . 4 

A 

V  2o 

. 60A02A3 , 1 • 117939?  1 • 1179398 

V  20 

1.1 <75398. 1.1 179398 

A 

A2 

A2 

iL 

S2T 

. 39AAAA ,.511 8055 ,.511 8g>5 
.5118055,  .5.-  180 A 5 

n 

A 

i, " 

.9062  , . 78AI66, . 78A166 

S2T 

.  78A  <  66  ?P’<166 

ES” 

. 020833,. 02Go3  ..020833 

ES2 

.C20833, .020833 

A 

C2L 

2A0.  ,2A0.  ,2A0.  ,2A0.  ,2A-J. 

50 


AFFDL-TR-71-155 
Part  1 1  I 


1 

12 

67 

INDC01 

1 

MSMAIN 

2 

CTABOI 

0., 2. 1,1. ,2. ,3. ,4. ,5. 

CTABOI 

3024. .3024. ,2160. ,2160.  ,2160.  ,2160. 

8 

CTABOI 

2160. ,2160. ,2160. ,2160. 

14 

I.N-DC02 

* 

CTAB02 

0.,  2. 1,1. ,2. ,4. ,5. 

CTAB02 

3024. ; 3024. ,2160. ,2160. ,2160. ,2160. 

8 

CTAB02 

2160. ,2160. ,2160. ,2160. 

14 

in:co3 

1 

NS2ND* 

2 

CTAB03 

0. 2. ,3. ,4  5. 

CTAB03 

72., 72., 72., 72., 72.  72., 72., 72., 72., 72. 

8 

1 NDCC4 

1 

CTABOI* 

0.,l.,i  ,2. ,3* ,4. ,5. 

CTABQ4 

720. ,7*0. ,720. ,720. ,720. 

8 

f„804 

,70. ,720. ,720. ,720. ,7*0 

13 

MASS  2 

3 . 1 08099 ,3.1 08099 ,3-1 08999 

i  »SS2 

3.108099,3.108099 

4 

MUS 

0. ,0. ,0. ,0. ,0.  , 

ES 

.02  333, .020833,. 020833 

ES 

.  020t.,3,  .02083: 

4 

SB 

1.8333. ,2. 08333, 2. 08333 

SB 

2.08333,2.08333 

4 

CASK 

1..0..0, 
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the  outside  tire  radius  rok  (Reference  1,  pg  110)  in  feet.  DELTAM  is 
an  array  of  maximum  tire  deflections  in  feet.  If  actual  tire  deflections 
exceed  these  values,  the  tire  is  assumed  to  blowout  and  the  program 
stops.  RLT  is  the  runway  length  in  feet.  IFD  is  an  indicator  which 
determines  whether  the  tire  force-deflection  curves  f ^ ( <Sk )  (Reference  1, 
pg  99)  are  input  as  tables  (i.e.,  table  FTAB01)  or  in  equation  form.  If 
I F0  has  the  value  one,  the  tire  force  deflection  curves  assume  the 
equation  form  a^  (6^)  b^.  The  ak  and  b^  constants  are  stored  in  the  Al 
and  Bl  arrays,  respectively.  The  units  of  the  data  in  the  Al  and  Bi  arrays 
must  be  such  that  when  <5^  is  in  feet,  the  function  is  in  pounds.  If  IFD 
has  the  value  zero,  the  tire  force  deflection  curves  are  contained  in  a 
single  two-dimensional  table  FTAB01.  This  table  setup  is  similar  to  .he 
thrust  table  TTAB10.  The  NSTRUT  indicator  (already  input)  is  synonymous 
with  the  I T 1  OX  indicator  of  TTAB10  and  indicates  the  number  of  force- 
deflection  curves  contained  in  FTAB01.  The  NDELTA  indicator  (not  shown) 
indicates  the  number  of  6^  points  in  each  cable  f^  (6^) .  The  NDELTA 
card  is  followed  by  the  FTAB01  table  which  lists  the  NDELTA  6^  points 
(in  ft)  In  ascending  order,  followed  by  the  gear  numbers  (i.e.,  1.,  2., 

3.,  etc.,  up  to  NSTRUT).  The  actual  force  data  (lbs)  for  gear  1  then 
follows  for  ascending  6^.  This  is  followed  by  the  force  data  for  gear  2, 
then  gear  3  etc.,  until  there  are  NSTRUT  tables  contained  in  FTABOI. 

FTAB02  is  the  runway  perturbation  table  £(XmJ  (Reference  I,  pg  110). 

The  independent  variable  is  X^  in  feet  and  the  dependent  variable  is 
the  runway  perturbation  c  in  feet  FTAB03  is  the  tire-runway  coefficient 
of  friction  table  (Reference  1,  pg  99)-  The  independent  variable  is 
"Pskdk"  ^e^erence  1 >  P9  218)  whose  range  is  0.  -  1.0.  All  tires  are 
assumed  to  have  the  same  U|<  properties.  Note  the  difference  in  format 
for  gear  array  input  and  gear  table  input  (pgs.  9'10). 

MOMENT  Is  an  array  of  the  moments  of  Inertial  1^  (slug  -  ft  )  of  a 
single  tire-wheel  combination  (Reference  I,  pgs  2 1 5- 2 1 6)  on  each  strut. 

MB  is  an  array  of  the  Initial  values  given  to  the  braking  moment 
MBk  ”  ft)-  The  values  in  the  MB  arroy  are  normally  zero  (Reference  1, 
pg  216).  The  RF  array  is  the  fully  extended  axle  position  rfk  (Reference  I, 
pg  82)  In  feet.  V 7  is  the  velocity  V0  (ft/sec)  In  Reference  I,  pg  135- 
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Its  value  should  be  z^ro.  The  PZERO  array  is  the  preload  pressure  PQk 
(lb/ft2)  of  the  upper  air  chamber  (Reference  1,  pg  86)  for  each  gear. 

The  VZERO  arrary  is  the  preload  volume  VQ ^  (ft3)  of  the  upper  air 
cnamber  (Reference  I,  pg  86)  for  each  gear.  The  A  array  is  the  area  A^ 
(ft2)  of  the  main  piston  (Reference  1,  pg  86)  for  each  strut.  The  P20, 
V20,  A2  arrays  are  the  preload  pressure  (ib/ft2),  preload  volume 

VQk2(ft^)»  and  piston  area  Aj^ft2)  (Reference  1,  pg  86),  respectively, 
'or  the  .  icondary  air  r'  umber  in  each  strut.  An  IL  indicator  value  of  0 
includes  the  scondary  oiston  and  air  chamber  in  the  gear-strut 
simulation.  An  IL  indicator  value  of  l  removes  the  secondary  piston  and 
a;r  chamber  from  the  jar-strut  simulation.  The  S2T  array  is  the  maximum 
allowed  displacement  (>t)  of  the  secondary  piston  (Referenct  1, 

pg  95“97)  for  each  strut.  The  ES2  array  is  the  integration  accuracy 
constraint  Eg,.,  (ft)  for  the  secondary  piston  position  of  each  strut. 

The  C2L  array  is  the  linear  drag  coefficient  Cj^l  ( lb-sec/ft)  of  the 
secondary  piston  (Reference  1,  pg  88)  for  each  strut.  The  CTAB01  table 
is  the  nonlinear  damping  coefficient  of  the  ma  in  orifice  C|<(lb-sec2/ft2) 
for  ma.n  strut  compression  (Reference  I,  pg  93)  for  each  strut.  The 
C  AB02  table  is  the  nun II near  damping  coefficient  of  the  main  orifice 
Ck( lb-sec2/f t2)  for  main  strut  extension  for  each  strut.  The  CTAB03 
table  is  the  nonlinear  damping  coefficient  of  the  secondary  orifice 
C|<2  ( lu*sec2/f t2)  for  main  strut  compression  (Reference  1,  pg  93)  for 
uch  strut.  The  CTAPO^  table  is  the  nonlinear  damping  coefficient  of 
the  secondary  orifice  Ck2  (  Ib-secVfi  )  for  main  strut  extension  for 
each  strut.  Each  of  the  CTAB  tables  is  two-dimensional  and  similar  in 
input  format  to  the  FTABOI  table.  The  two  independent  variables  are 
strut  position  and  strut  number.  The  CTAB  tables  each  begin  with  the 
strut  deflections  followed  by  the  strut  numbers.  The  remainder  of  the 
data  in  each  table  Is  the  value  of  the  orifice  coefficient  for  increasing 
strut  deflection  beginning  with  strut  I.  NSMAIN  is  the  number  of 
independent  strut  positions  stored  for  each  main  strut  in  tables 
£TAB0I  and  CTAB02.  NS2NDY  Is  the  number  of  Independent  secondary  piston 
positions  stored  for  each  strut  In  tables  CTAB03  and  CTABOA.  If  a  table 
is  to  be  used,  its  associated  indicator  (i.e.,  INDC01,  etc.)  Is  needed. 
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The  MASS2  array  is  the  mass  (slugs)  of  each  secondary  piston 
(Reference  1,  pg  92)  in  each  strut.  The  array  MUS  is  the  coefficient  of 
sliding  friction  at  the  wing  gear  root  (i.e.,  where  the  main  strut 

is  supported  by  the  wing)  for  each  strut  (Reference  1,  pg  88).  The  ES 
array  is  the  integration  accuracy  constraint  E^  (ft)  for  the  main  strut 
position  (Reference  1,  pg  95"97)  of  each  strut.  The  SB  array  is  the 
maximum  allowed  displacement  (ft)  of  the  main  strut  (Reference  l, 
pg  95-97)  for  each  strut.  The  CASK  array  is  an  array  of  indicators 
denoting  whether  or  not  the  corresponding  gear  is  castered  (i.e.,  its 
wheels  allowed  to  align  themselves  with  the  axle  velocity  vector). 

A  one  value  causes  the  gear  to  be  castered.  A  zero  value  fixes  the 
orientation  of  the  axle  relative  to  the  body  axes  system. 

3.  AUTOPILOT  DATA 

Column:  1  12 

I NDAUT  1 


A  zero  value  for  the  INDAUT  indicator  deletes  all  autopilot 
calculations;  a  one  value  causes  the  autopilot  calculations  to  be  made. 


a.  Engine  Data 


Column:  1  12 


IN  * 

ZN  5. *1.3. 39. 3. 39, 5. 41 

YN  61.8,39.7,-39.7,-61.8 

N  I .975.-1 .. -I .,1 .975 


The  IN  indicator  Is  used  to  indicate  the  number  of  engines  on  the  aircraft 
(Reference  1,  pg  197).  A  minimum  of  one  engine  and  a  maximum  of  four 
engines  can  be  simulated.  The  ZN  and  YN  arrays  are  the  Z,  Y  body  coordi¬ 
nates  (In  ft)  of  each  engine  thrust  vector  origin.  Because  of  the  throttle 
autopilot,  the  data  In  the  ZN  and  YN  arrays  must  correspond  to  the  engine 
physical  arrangement  In  one  direction  along  the  wing.  As  discussed  In 
Reference  I  (pg  227),  the  thrust  vectors  are  assumed  parallel  to  the 
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X  longitudinal  body  axis;  therefore,  the  engine  X  coordinates  are  not 
needed.  The  N  array  contains  the  initial  throttle  settings  for  each 
engine. 


b.  Drag  Chute  Data 


Column:  1 


ICS  the  drag  chute  indicator  signal  (Reference  1,  pgs  63  6  179),  is 
usually  initialized  at  a  zero  value.  The  variables  CDCH,  SSH,  XCH,  YCH 
and  ZCH  are  the  variables  C^^,  (ft^)  ,  X^  (ft),  YCH  (ft)  and 

ZCH  (ft)  respectively  (Reference  1,  pgs  61  S  62). 


c.  Phase  Begin  -  Terminate  Data 


Column:  1  12 


ITO 

0 

HF 

100. 

NF 

0 

XRF 

100. 

HRF 

25. 

DELTAH 

3. 

KP 

0 

T 1 

0. 

NLRI 

1 

VS 

10. 

XS 

10000. 

TS 

40. 

HS 

0. 

The  I  TO  Indicator  Is  the  takeof f- land! ng  iidicator  (Reference  1, 
pgs  180  t  182).  A  zero  value  Is  used  i or  a  landing  problem  and  a  one 
value  for  a  takeoff  problem,  HF  is  the  flare  altitude  hp  (ft)  (Reference 
pg  182).  NF  Is  the  HF  glide  slope  termination  indicator  (see  Reference  1 
pg  182).  An  NF  value  of  zero  allows  the  problem  to  go  into  the  flare 
phase  when  the  hp  altitude  Is  reached.  An  Np  value  of  one  terminates 
the  problem  at  the  glide  slope  end  (l.e.,  hr  altitude).  The  variables 
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XRF  and  HRF  are  the  initial  values  of  the  variables  Xgp  Cft)  and  hRp  (ft) 
(Reference  1,  pgs  164  and  182).  They  are  primarily  used  to  start  the 
problem  in  the  "hold  mode".  DELTAH  Is  the  variable  6^  (ft)  (Reference  1, 
pgs  180  £  182).  Because  of  possible  high  flare  guidance  acceleration 
(i.e.,  A^g)  as  flare  termination  is  reached,  6^  should  always  have  a 
nonzero  positive  value.  KP  is  an  indicator  (Reference  1,  pgs  1 80  £  l 82) 
used  to  start  the  problem  in  the  landing  roll  phase.  A  KP  value  of  zero 
allows  the  program  to  establish  its  own  time  of  impact,  T|,  and  a  KP 
value  of  one  requires  the  impact  time  T|  (sec)  to  be  read  into  the 
program  as  Tl.  NLRI  is  the  NLR  Indicator  (Reference  1,  pgs  180  f.  182). 

An  NLRI  value  of  one  stops  the  problem  on  impact,  and  a  value  of  zero 
allows  the  problem  to  continue  after  impact  and  terminate  on  one  of  the 
three  rurway  stopping  condition  Vg  (ft/sec),  X5  (ft),  T$  (sec)  (Reference  1, 
pg  182)  which  are  input  as  VS,  XS,  and  TS,  respectively.  HS  is  the 
runway  altitude  hs  (ft),  at  which  the  takeoff  roll  is  terminated 
(Reference  I,  pgs  178  £  181). 


d.  Takeoff  Condition  Data 


Column : _ 1 _ 1_2 

VATO  0 . 

ALPHT }  0. 

VATO  and  ALPTHO  are  takeoff  conditions  V^q  (ft/sec)  and  Oj-q  (deg) 
respectively  (Reference  1,  pgs  1 78  and  1 8 1 ) . 

e.  Glide  Slope  Maneuver  Data 


Col  umn : _ 1  1 2 

VD  220. 

DELVE  0. 

EPSGS  3. 

ALPHOS  -5. 

ALPHDL  M*. 

HCG  20.56 

DELEPS  .05 

RFH  5. 

PGS  .04 

DELS IG  .1 

f;FY  10. 

PH  1C  1.0 
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VO  if  the  desired  inertial  velocity  down  the  glide  slope  V^  (ft/sec) 
(Reference  1,  pg  152).  DELVE  is  the  velocity  error  A  Ve  (ft/sec)  which 
shojld  be  input  as  zero.  EPSGS  is  the  glide  slope  elevation  (deg) 
(Reference  1,  pg  148).  ALPHDS  .  nd  ALPHDL  are  the  lower  and  upper  limits 
afiS  (deg)  and  adL  (deg),  respectively  (Reference  1,  pg  160)  for  the  glide 
slope  desired  angle  of  attack.  HCG  is  the  variable  h^Q  (ft)  (Reference  1, 
og  148).  DELEPS  is  the  allowed  vertical  glide  slope  angular  error  6£  (deg). 
RFH  is  the  rate  feedback  constant  RF^  (sec)  and  PGS  is  the  phugoid  control 
constant  PGS  (deg/ft)  (Refe> ence  1,  pgs  151  &  157).  DELSIG  is  the  allowed 
horizontal  glide  slope  anguiar  ;rror  <5a(deg),  RFY  is  the  rate  feedback 
constant  RF,(sec),  and  PH  1C  is  the  roll  control  angle  4>c  (deg) 

(Reference  1,  pgs  149,  15',  l Svi) . 

f.  Flare  Maneuver  Oc ta 


1 

12 

XTD 

ICO. 

HTD 

25- 

VXTD 

200. 

VHTD 

-1  , 

1  IR 

0 

LD 

5000. 

DA 

.02 

TL 

5000. 

TU 

126000. 

XTD  is  the  variable  (ft)  which  locates  the  runway  beginning  with 
respect  to  the  glide  slope  origin  (Reference  1,  pg  167).  HTD,  VXTD  and 
VHTD  are  the  In  tial  desired  touchdown  conditions  h-jp  (ft),  V^p  (ft/sec), 
and  VhTD  (ft/sec I,  respectively.  MR  is  the  Indicator  IR  (Reference  1, 
pg  173)  and  shoi  Id  always  have  a  zero  Initial  value.  LD  is  the  expected 
runway  landing  distance  lQ  (ft).  DA  is  the  accuracy  (in  deg)  of  the 
search  done  t<>  solve  the  nonlinear  simultaneous  equations  (Reference  1, 
pgs  171  6  176,  Equations  299  and  300)  for  the  flare  ad ,  Td  comands. 

TL  (lbs)  and  TU  (lbs)  are  the  lower  and  uppe  bounds,  respectively,  on 
the  desired  thrust,  T^,  in  the  flare. 
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g.  Hold  Maneuver  Data 
Column:  1  12 


ALPDES 

TTD 

KF 

PM 


6.5896 
*♦7073. 
INT  1,1 ,1,1 
1  1  . 


If  the  problem  is  started  in  the  hold  mode,  (deg)  and  Td  (lbs) 
must  be  input  to  the  problem  as  the  variables  ALPDES  and  TTD, 
respect i valy.  KE  is  the  "kill  engine"  indicator  array  (Reference  1, 
pgs  172,  I98,  &  203).  This  allows  the  option  to  pull  the  throttles  back 
to  forward  idle  (if  the  engines  are  not  in  a  fixed  mode)  in  the  hold  mode. 
A  one  value  in  the  array  exercises  the  kill  power  option  for  that 
particular  engine,  and  a  zero  value  leaves  that  particular  engine  at  the 
previous  fixed  thrust  level.  PM  is  the  tail  down  constraint  angle  Pm  (deg) 
(Reference  1,  pgs  172  6  177). 

h.  Landing  F;ol  1  Maneuver  Data 


Column:  1 


12 


TSP 

TRV 

TCH 

TBK 

ISS 

ILR 

IBS 


1  . 
3. 
0. 
2. 
0. 
0 
0 


TSP,  TRV,  TCH  and  TBK  are  the  variables  tQn  (sec),  t  (sec),  t 


sp 


ch 


(sec),  and  t..  (sec),  respectively  (Reference  1,  pg  179).  ISS,  ILR,  and 

Dk 

IBS  are  the  initial  values  of  the  indicators  ISS,  ILR,  and  IBS,  respec¬ 
tively  (Reference  1,  pg  179).  These  indicators  are  normally  input  as  zero 
(note  ISS  is  the  only  indicator  that  is  floating  point). 
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Engine  Failure  Stage  Data 


Column:  1  8  12 


1C 

INT 

1,1 

JJ 

XRF1 

0. 

IT1 

INT 

1  ,  1 

.1,1 

XRF2 

0. 

IT2 

INT 

1,1 

,1,1 

HI 

0. 

1  H 1 

INT 

1 , 1 

,1,1 

H2 

0. 

IH2 

INT 

1 , 1 

,1,1 

HR1 

0. 

IHR1 

INT 

1 . 1 

,1,1 

HR2 

0. 

IHR2 

INT 

1 , 1 

,1,1 

TR1 

0. 

ITRI 

INT 

1 , 1 

,1,1 

TR2 

0. 

ITR2 

INT 

1 , 1 

,1,1 

The  1C  array  contains  the 

initial  values  for  the  1C  (l)  array 

(Reference  1 ,  pgs 

183  6 

184)  . 

A  one  value  in  the  array  indicates  the 

engine  is  working, 

and  a 

zero 

value  indicates  engine  failure.  Reference 

pg  184  and  the  following  indicates  the  other  variables: 


XRF 1 

-  xRF, (ft) 

IT1 

-  IT  1(1) 

XRF2 

xRF2(ft) 

IT2 

-  1 T2  ( 1 ) 

HI 

-  h 1 (ft) 

1  HI 

-  IH1(I) 

H2 

-  h  ^ (ft) 

IH2 

-  182(0 

HRI 

-  hR)  (ft) 

1  HR  1 

-  [HR  1 { 1 ) 

HR2 

-  hp2(ft) 

IHR2 

-  1 HR2  ( 1  ) 
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TR 1  -  t  j (sec) 

ITR1  -  ITRl(l) 
TR2  -  tr2(sec) 
ITR2  -  ITR2 ( I ) 


j.  Brake  Condition  Stage  Data 


Column:  1 

8 

12 

IB 

INT 

0,1, 1,1,1 

TBK1 

100. 

IBKI 

INT 

0,1, 1,1.1 

TBK2 

100. 

IBK2 

INT 

0,1, 1,1,1 

The  IB  array  contains  the  initial  values  of  the  brake  condition  array 
Ig.  (Reference  1,  pgs  183  £  185).  The  variables  TBK1  and  TBK2  are  t^ 
(sec)  and  tbk2  (sec),  respectively.  The  arrays  IBKI  and  IBK2  are  the 
arrays  *Bkl(l)  and  lBk2(l),  respectively. 


k.  Pitch  Autopilot  Data 


Column:  1 

12 

1 

12 

ALPDL 

1  . 

DFLQC2 

0. 

RFALPH 

.5 

TST 

5.0 

DELALA 

.01 

DELQF 

0. 

PSH 

-20. 

DELFD1 

1 . 

PSH2 

-20. 

DELQT0 

0. 

RFALP2 

.5 

DELQL 

-15 

DFLQC 

0. 

DELQll 

25. 

ALPDL  (deg/sec)  is  the  maximum  ad  allowed  in  the  pitch  autopilot 

(Reference  1,  pg  191)  RFALPH  is  the  rate  feedback  constant  Rp  (sec) 

(Reference  I,  pgs  188  £  189).  DELALA  is  .the  allowed  angular  error 

Aaa  (deg).  PSH  is  the  overcontrol  sensitivity  PS^  (deg/deg)  (note  sign). 

RFALP2  and  °SH2  are  the  flare,  hold,  landing  roil  C  takeoff  values  for 

RF  and  PSU,  respectively  (Reference  1,  pg  189).  DELQC  and  DELQC2  are 
a  h 

the  "bang-bang"  control  magnitude  6  (deg)  (Reference  I,  pg  189  6  191) 
for  glide  slope  and  flare  (etc),  respectively.  TST  is  the  nose-over 
time  t  (sec)  (Reference  I,  pgs  188  l8y).  DtLFDI  is  the  angular 
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rate  6p  (deg/sec)  and  DELQF  is  the  final  nose-over  elevator  deflection 

0qp  (deg)  (Reference  1,  pgs  1 88  £  189).  DELQTO  is  the  takeoff  elevator 

initial  deflection  (deg).  DELQL  and  DELQU  are  the  lower  anj  upper 

elevator  deflection  limits  6  .  (deg)  and  6  (deg),  respectively 

C|  U  C|  u 

(Reference  1,  pgs  1 88  £  190). 


1.  Yaw  Autopilot  Data 


Column:  1 _ 1_2 

RFB  .5 

DELBA  .05 

PSR  10. 

RFPSI  .5 

DPS  I A  .05 

PSPSI  -10. 

DELRL  -35- 

DELRU  35. 


RFB  is  the  rate  feedback  constant  Rpg(sec)  (Reference  1,  pgs  192  £  193) 
DELBA  is  the  allowed  angular  error  Ag  (deg).  PSR  is  the  overcontrol 

d 

constant  PSR  (deg/deg).  RFPSI  is  the  rate  feedback  constant  Rp^  (sec) 
(Reference  1,  pgs  192  £  193).  DPSIA  is  the  allowed  angular  error  Aip  (deg). 

3 

PSPSI  is  the  overcontrol  constant  (note  sign)  PS^  (deg/deg).  DELRL  and 
DELRU  are  the  lower  6rL  (deg)  and  upper  <$ru  (deg)  rudder  deflection  limits 
respect i vel y . 

m.  Roll  Autopilot  Data 


Column:  1 _ 1_2 

RFPHI  .5 

DPHIA  .05 

PSA  -15. 

DELPL  -60. 

DELPU  60. 


RFPHI  is  the  rate  feedback  constant  Rp^  (sec)  (see  Reference  1, 
pgs  191*  £  195).  DPHIA  is  the  allowed  angular  error  A<J>  (deg) .  PSA  is  the 

3 

control  constant  PS^  (deg/deg).  DELPL  and  DELPU  are  the  lower  6p^  (deg) 
and  upper  6p^  (deg)  aileron  deflection  limits  respectively. 
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n.  Throttle  Autopilot  Data 


Co  1 umn : 

1 

8 

12 

67 

TF 

INT 

0,1, 1,0 

NDF 

0. ,-i . ,-l . ,0. 

IR 

INT 

1,1, 1,1 

NB 

1 .05,1.05,1.05,1.05 

NLR 

-2.  ,-2. ,-2. ,-2. 

NTO 

2. ,2. ,2. ,2. 

K2 

.5 

K3 

•  5, .5, .5, .3333, .3333 

K** 

•  5,.  5,.  5,. 5,. 3333,. 3333,- 3333,. 3333 

K** 

.25, -25, .25 

9 

The 

TF  array 

is  the 

throttle  fix  indicator  arrcy  TF  (l) 

(Reference  1  , 

pg  20**).  The  NDF  array  is  the  fixed  throttle  setting  array  N^p.(l).  The 
IR  array  is  the  engine  reverse  capability  indicator  array  IR  (I) 


(Reference  1,  pgs  198,  20**,  208,  210,  212-215).  The  NB  array  is  the  engine 
throttle  constraint  array  for  reverse  Ng  (l).  The  NLR  array  is  the  reverse 
throttle  setting  array  for  landing,  N^  (I).  The  NTQ  array  is  the  takeoff 
throttle  setting  array  N^.g  (l).  K2  is  the  engine  load  factor  ^ ( 2 )  1 2 1 
(Reference  1,  pgs  203  6  20**).  The  K3  array  contains  the  engine  load 

factor  constants  k(3)131>  k(3) 232 '  k (3) 121 »  k(3) 1231 '  k(3)1232’  1 n  that 
order.  The  K**  array  contains  the  engine  load  factor  constants 

k  (*♦)  232  ’  k(**)3**3’  k(**)  2**2  ’  k  ( *4 )  2  3*+2  ’  k(**)23**3’  k(**))3**l’  k  (**)  1 3**3  ’ 
k(**)  123**'  ’  k (**)  1 23**2 ’  k(**)l23**3’  in  that  order> 


o.  Brake  Autopilot  Data 


Column : 


1 


12 


HBC  0. ,0. ,0. ,0. ,0. 

PD  .01 

DELTAW  .001 

0MECDI  20. 

MBL  0. ,0. ,0. ,0.,0. 

MBU  3.E5,3.E5,3.E5,3.E^,3.F.5 


The  MBC  array  is  the  constant  braking  array  Hgr|*lb-ft)  (Reference  1, 


pgs  217  6  218.  PD  is  the  deslr. 


percent  stod  P^.  DELTAW  is  the  allowed 


fraction  wheel  speed  error,  Ao.  <j)MECDl  is  the  control  tire  angular 
acceleration  w.  rad/ _ 2  •  The  arrays  MBL  and  PBU  are  lower  nlgL.  (lb-ft) 


and  upper  H 
21°  r,  222) 


Bu  i 


c  sec 

(lb-ft)  braking  moment  limit  arrays  (Reference 


P9S 
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p.  Control  Response 

Data 

Column:  3 

12 

DELHS 

25. 

DELRRD 

35. 

DELA 

60. 

NED1 

.125 

OELHS,  DELRRD,  DELA,  and  NED  1  are  the  control  variable  rates 

^(^sec)’  ^DV^sec)*  6A  (^sec)’  and  ^(^sec)’  respectively 
(Reference  pg  224)  . 

q .  Ini t i al i zat ion 


Column:  1 

12 

1 

12 

IAP 

1 

MANLOG 

1 

HR 

0. 

PITCHP 

1 

DELQN 

0. 

YAWAUP 

1 

DELQDE 

.07227 

ROLLAP 

1 

DELQD 

.07227 

THROAP 

1 

DELPD 

0. 

BRAKAP 

1 

DELRD 

0. 

CONTRP 

1 

AUXICP 

1 

INDICP 

1 

IAP,  HR,  DELQN ,  DELQDE ,  DELQD,  DELPD ,  and  DEI  RD  are  initial  input 
values  of  IAP,  hR(ft),  Sqn(deg),  6qd(deg),  6q(deg),  6p(deg),  6r(deg)  , 
respectively  (Reference  1,  list  of  symbols). 

r.  Ac.opilot  Output  Indicators 

There  are  nine  categories  of  autopilot  output:  auxiliary  computations, 
maneuver  logic,  pitch  autopilot,  yaw  autopilot,  roll  autopilot,  throttle 
autopilot,  brake  autopilot,  .ontrol  response,  and  status  indicators.  The 
following  nine  indicators:  AUXICP,  MANLOG,  PITCHP,  YAWAUP ,  ROLLAP, 

THROAP,  BRAKAP,  CONTRP,  INDICP  are  associated  iespectively  to  the  nine 
categories  of  autopilot  output.  The  Indicators  initially  have  a  zero 
value  which  eliminates  all  autopilot  output.  If  the  indicator  on  input 
is  given  a  one  value,  the  output  wilt  be  printed.  See  Section  III, 
paragraph  2,  for  details.  If  used,  these  indicators  follow  the  initial¬ 
ization  card  inputs  just  discussed. 
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s.  Plot  Tape  Data 

A  magnetic  ‘■ape  may  be  used  to  output  selected  rigid  uody  data  and 
selected  landing  gear  data.  A  request  or  label  card  is  needed  after  the 
J4>B  card  to  mount  a  tape.  The  following  data  is  cJso  required  in  the 
first  stage  of  the  data  decks. 


Card  Column:  1 _ 1  2 

I PLT  I 
ISDF  1 
ISTPL1  1 
ISTPL2  1 
ISTPL3  1 
ISTPU  1 
ISTPL5  1 


IPLT  =  1  denotes  that  data  will  be  saved  on  a  tape  for  plotting. 

ISDF  =  1  denotes  rigid  bod«  data  will  ‘e  saved  on  tape.  1STPL1  *  l, 
ISTPL2  =  1,  ISTPL3  =  1,  ...  etc.,  denote  that  data  for  landing  gears 
numbers  1,  2,  3,  ...  etc.,  will  be  saved  on  tape.  ISDF  =  1  saves  Lm, 
Mm»  ^m*  q,  6p,  typ,  4>p»  Fzm*  3.,*  *g»  hg,  Xg,  Yg,  ax,  Xg  at  each  printed 
time.  ISTPLk 
S„,  ,  M 

i 

needed  to  further  reduce  this  data  for  Calcomp  plot.  .  Use  of  this 
separate  program  is  explained  in  the  last  volume  of  this  repo't  . 


1  saves  FTk,  Spk,  $k,  P^.  ^k’  $k .  $k>  $2k,  $2k, 

2k’  "Ak’  uTk '  ^k  at  each  Panted  time.  A  separate  computer  program  is 


A.  STAGING  DATA 

Pages  32~33  explain  the  basic  staging  options  available.  One 
additional  staging  option  was  added  whicn  will  be  explained  later.  The 
general  staging  logic  is  so  built  that  almost  any  kind  of  staging  can  be 
done.  As  shall  be  seen,  the  specific  staging  done  in  ’arge  part  controls 
the  time  efficiency  of  the  program  and  parts  of  the  staging  is  a  must  to 
even  get  th»  program  to  run. 


a.  Staging  Gears  Into  Program 


The  added  calculations  to  account  for  the  landing  gears  is  voluminous 
and  time  consuming.  Even  though  the  gea  computations  can  be  made  and 
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correct  answers  obtained  during  \  tide  free  flight  (i.e.  ie  slow 
and  fiare),  these  calculations  are  unnecessary.  The  following  staging 
is  therefore  done  first: 


Column : 


1 

8 

12 

INDLG 

0 

1  STAGE 

0 

DECRES 

BCD 

1  HR 

STESTD 

TRA 

?5- 

INDLG 

-1 

particular  value  of  hR  cf>r  this  s  ge  must  be 


The  zero  INOLG  indication  prevents  the  geo  •  calculations  from  being 

made.  As  the  aircraft  nears  the  ground  and  impact  approaches  (this  is 

sensed  by  a  decreasing  test  on  the  runwU(  altitude,  (ft),  the  gear 

calculations  are  staged  into  the  program.  In  this  particular  tr  t,  the 

gear  calculations  were  su.ged  in  whe  h^  <  25  ft.  Th  hange  in  data 

required  is  a  value  f  -1  for  the  INDLG  indicator.  Note  that  the 

fficiently  greater  than 

h  to  insure  that:  impact  does  not  occu-  before  the  stane  is  made.  When 

eg 

the  gears  are  staged  in,  the  init'  '  values  of  main  str^k  tosition,  S^, 
main  strut  velc  ity,  S^,  s  ondary  piston  position,  Sk2’  secondary 
piston  velocity,  S^,  and  tire  angular  rate,  are  all  automatically 

zero.  If  the  problem  begin:  on  the  runway,  as  in  the  takeoff  roll,  the 
initial  values  in  the  arrays  3^ ,  S^,  5^,  S^2  and  ui^  must  be  read  in. 
This  is  done  through  tie  following  variables  which  are  placed  immediately 
after  the  INDLG 


!  card: 

SI 

S, 

(ft) 

S  -i  2 

S2 

(f  t ' 

S3 

S3 

(ft) 

sk  ~ 

S<, 

(ft) 

S5 

S5 

(ft) 

SD1 1  ~ 

S, 

(ft/sec) 

bD12  ~ 

S2 

(ft/sec) 
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SDI 3 

S3 

(ft/sec) 

SD1 A 

sw 

S4 

(ft/sec) 

SDI  5 

/V 

S5 

(ft/sec) 

S21 

S12 

(ft) 

S22 

S22 

(ft) 

S23 

/%/ 

S32 

(ft) 

S24 

S^*2 

(ft) 

S25 

<v 

S52 

(ft) 

S2D1  1 

*12 

(ft/ sec) 

S2D12 

S22 

(ft/sec) 

S2013 

*32 

(ft/sec) 

S2D1 4 

S42 

(ft/sec) 

S2D15 

»v 

*52 

(ft/sec) 

0MT1 

WT1 

(rad/sec) 

0MT2 

(rad/sec) 

0MT3 

/N* 

UT3 

(rad/sec) 

OMTA 

“t«. 

( rad/sec) 

0MT5 

rw 

^5 

(rad/sec) 
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This  gear  stage  is  always  the  first  stage  and  is  considered  a  must 
for  a  time  efficient  program. 

b.  Smooth  Impact  Stage 

The  second  stage  (also  a  must)  can  be  looked  up  as  an  impending 
impact  stage. 


Column:  1 

8 

12 

A1NCRS 

BCD 

LDDELT10DELT2DDELT60DELT5 

STEST 

TRA 

-l.,-l.,-l.,-1. 

PRINT 

.01 

DELTS 

.01 

AMAXER 

.001 

PRTMIN 

0. 

The  heights  of  each  wheel  bottom  surface  (see  <5^  'n  Appendix  It) 
above  the  runway  (note  this  is  a  negative  value)  are  given  the  names 

DDELTl  -  DDELT5.  As  these  numbers  approach  zero,  impact  occurs.  On 

this  particular  test,  the  stage  was  performed  when  bottom  surface  of 
wheels  1,  2,  4,  or  5  was  less  than  one  foot  off  the  runway.  At  this 
point  the  required  data  change  is  a  small  VSRK  maximum  interval,  that 
is  AMAXER.  This  insures  a  smooth  impact.  Note  that  when  AMAXER  is  changed, 
the  integration  data  constraint  PRTMIN  5  AMAXER  S  DELTS  £  PRINT 
must  still  hold.  Note  also  that  the  h^  value  in  the  initial  stage  must 
be  large  enough  to  insure  that  the  wheel  bottom  surfaces  are  at  least  one 
foot  off  the  runway  (in  this  case)  when  the  gears  are  stage  in.  The 

reason  for  this  is  that  the  variables  DDELTl - DDELT5  don't,  even  exist 

until  the  gear  calculations  are  staged  in.  This  is  a  good  time  to  assert 
an  important  point  about  the  staging  logic:  that  is,  the  stages  will 
only  occur  In  the  order  they  appear  in  the  data  deck.  Therefore,  one 
must  make  sure  that  the  stages  are  In  the  proper  sequence  and  staged 
on  the  variables  appropriate  to  sense  that  sequence, 

c.  Efficient  AMAXER  Stage 

t 

The  small  AM^AXER  of  the  previous  stage  wa-.  used  solely  to  insure  a 
smooth  impact,  To  leave  AMAXER  at  this  small  value  would  not  be  very 

N  T 
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time  efficient,  thereby  severely  limiting  the  capability  of  the  VSRK 
integration  interval-  As  soon  as  impact  occurs  occurs  and  the  VSRK 
technique  has  had  an  opportunity  to  pick  an  appropriate  integration 
interval,  AMAXER  should  be  staged  back  to  a  reasonable  value.  This  is 
done  as  follows: 

Column:  1  8  12 


AINCRS 

BCD 

4DELTA J  DELTA2DELTA4DELTA5 

STEST 

TRA 

PRINT 

.08 

AMAXER 

,0k 

DELTS 

.0A 

The  tire  deflections  OELTAI  -  DELTA5  are  used  as  the  staging 

variable.  The  particular  tire  deflection  in  this  stage  was  0.1  foot. 

d.  Spoiler  Aero  Stage 

If  the  aircraft  is  landing,  the  next  stage  is  spoiler  actuation 
(if  spoilers  exist  on  the  aircraft).  The  staging  variatle  is  the  spoiler 
actuation  signal  ISS.  The  required  data  change  is  the  appropriate 

Column:  1  8  12 


AINCRS 

BCD 

1  ISS 

STEST 

TRA 

1. 

ATAB10 

o 

r*"k 

O 

O 

ro 

O 

ATABI 1 

.1162, .1162 

ATABI2 

0.  ,0. 

ATAB51 

.3^39.-3^39 

ATAB52 

-.0327. -.0327 

ATAB53 

0.  ,0. 

ATAB80 

.267.-267 

aerodynamic  coefficient  tables  with  spoilers  activated.  Since  the  spoiler 
signal  is  based  on  a  specific  time  after  timpact,  t  ,  one  has  to  make 
sure  that  the  staging  tire  deflection  of  the  previous  stage  is  not  too 
large  so  as  not  to  occur  before  t  .  If  tjp  is  very  small  (l.e.,  like  0. 
etc.),  then  the  spoiler  stage  may  have  to  be  placed  before  the  efficient 
AHAXER  stage. 
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e.  Bounce  Stage 

A  combination  of  high  sink  rate,  no  spoilers  (or  late  spoiler 
actuation),  stiff  landing  gear  strut,  and  a  light  aircraft  can  create  a 
situation  where  the  entire  aircraft  may  bounce  completely  off  the 
runway.  To  sense  this  condition,  the  staging  logic  had  to  be  modified. 
By  giving  the  ISTAGE  indicator  a  value  of  one  (it  was  initially  zero), 
all  four  variables  (not  just  one)  in  the  stage  testing  logic  must  meet 
the  increasing  or  decreasing  constraints  before  the  stage  is  performed. 
The  bounce  stage  is  therefore  as  follows: 

Column:  1  8  12 


1  STAGE 

1 

DECRES 

BCD 

4DELTA1 DELTA2DELTAkDELTA5 

STESTD 

TRA 

0. ,0. ,0. ,0. 

AMAXER 

.001 

DELTS 

.01 

PRINT 

.01 

1  STAGE 

0 

AINCRS 

BCD 

kDDELTlDDELT20DELTl4DDELT5 

STFST 

TRA 

.1..1..1..1 

AMAXER 

.05 

DELTS 

.  1 

PRINT 

.  1 

First  a  decreasing  test  is  done  to  check  if  the  aircraft  bounced; 
that  is,  see  if  the  tire  deflections  DELTA1  —  DELTAS  have  reached 
zero.  AMAXER  is  then  given  a  smai 1  value,  as  before,  to  Insure  a  second 
smooth  impact.  ISTAGE  is  given  a  zero  value  again  and  a  second  stage 
similar  to  the  efficient  AMAXER  stage  follows.  This  staging  Is  sufficient 
to  get  the  aircraft  through  one  complete  bounce.  As  many  bounce  stayes 
as  one  feels  Is  needed  can  be  used. 

f.  Data  Deck  End 

The  last  stage  is  ended  with  the  following  cards: 

Column  : _ 1  8 _ 

INDSTF  I 

TRA 

END  of  JOB  card 

This  concludes  the  data  deck  setup. 
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SECTION  VI 
OUTPUT 

The  time  intervals  for  computer  output  are  PRTMIN  and  PRINT  (see 
discussion  on  pages  47-48).  As  previously  discussed,  the  PRTMIN  option 
should  not  be  used  unless  one  is  really  interested  in  the  high  frequency 
response  of  the  landing  gears.  The  long  term  response  of  the  gears  and 
main  airframe  can  be  approximated  by  a  PRINT  value  of  about  0.1  or  0.2 
seconds.  The  selection  of  a  PRINT  value  should,  however,  be  tempered 
by  the  frequency  of  the  specific  output  to  be  examined.  The  reason  for 
this  print  interval  caution  Is  that  there  can  be  as  much  as  50  lines  of 
output  each  time  a  PRTMIN  or  PRINT  option  is  exercised.  For  the  computer 
system  available  to  AFFOL,  experience  has  shown  the  TOLA  problem  to  be 
output  bound  In  such  a  case.  Therefore,  use  of  a  PRINT  interval  smaller 
than  that  which  is  needed  will  result  In  a  computation  time  longer  than 
needed. 

1.  TRAJECTORY  PRINTING  METHOD 

The  printing  of  a  trajectory  may  be  divided  into  four  categories. 

a.  Initial  Printing  -  The  printing  of  specific  values  at  the  first 
stage  and  at  each  subsequent  major  stage. 

b.  Code  Printing  -  The  printing  of  codes  which  will  identify  the 
variables  which  are  to  be  obtained  in  the  coming  time  history  print. 

c.  Time  History  Printing  -  The  printing  of  values  specified  at  the 
requested  points  of  the  trajectory. 

d.  Diagnostic  Error  Printing  -  The  printing  of  errors  oetected  by 
the  program. 

All  input  data  involved  for  a  case  is  printed  on  the  output  page 
preceding  the  computation  of  the  first  stage  printout.  Also,  data  read 
in  at  stage  times  will  be  printed  out  between  the  stages  of  the  trajectoiy 
output . 
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Initial  print  is  designed  to  print  certain  value  which  will  be 
constant  during  the  trajectory  and  series  as  a  reminder  of  what  values 
have  been  uaed  for  these  constants. 

Code  printing  is  performed  once  per  major  stage  o  identify  the 
time  h is tory . 

The  time  history  print  is  designed  to  print  in  a  minimum  space. 

That  is,  if  a  certain  variable  is  not  desired  as  output,  it  is  not 
printed  and  other  desired  variables  are  moved  in  the  print  format 
accordingly. 

The  entire  printing  is  controlled  to  print  on  a  page  11  x  14  inches 
and  will  print  a  maximum  of  5^  lines  per  page.  Page  ejection  and  line 
control  are  provided  by  the  subroutines  DEF  and  LINES. 

There  are  three  basic  kinds  of  TOLA  output: 
main  airframe,  autopilot,  and  landing  gears. 
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2.  MAIN  AIRFRAME  OUTPUT 

There  are  two  places  during  the  run  which  may  print  a  set  of  symbols, 
initial  print  and  time  history  print.  This  section  will  present  (1)  the 
computer  symbol,  (2)  its  associated  engineering  symbol  (if  any), 

(3)  the  unit  of  measure,  and  (A)  remarks. 

( 1 )  Initial  Print 

None 


(2)  Time  History  Print 


Symbol 

Quant i ty 

Uni  ts 

Remarks 

ALPHD 

a 

deg 

Angle  of  attack 

ALPHD1 

a 

deg/sec 

Rate  of  change  of  angle  of  attack 

AMACH 

mn 

Mach  Number 

AX77F 

ax 

ft/sec2 

Body  component  of  inertial  acceleration 

AXP7F 

Axp 

ft/sec2 

Platform  -»cc>ieration 

AY77F 

ay 

ft/sec2 

Body  co*'  nent  of  inertial  acceleration 

AYP7F 

Ayp 

ft/sec2 

Piatfc  acceleration 

AZ77P 

az 

f t/sec2 

Body  component  of  inertial  acceleration 

AZP7F 

Azp 

f t/sec2 

Plat  am  acceleration 

BETAD 

8 

eg 

Angl i  )f  s ides  1 i p 

BETAD1 

8 

deg/sec 

Rate  of  change  of  ang^e  of  sideslip 

DYNPP 

q* 

lbs/ft2 

Dynamic  pressure 

FCX 

Fc* 

1  bs 

Body  component  of  drag  chute  force 

FCY 

Fc/ 

lbs 

Body  component  of  drag  chute  force 

FCZ 

1  bs 

Body  component  of  drag  chute  force 

FDC 

fdc 

lbs 

Total  drag  chute  force 

GAM  70 

Y 

deg 

Elevation  flight  path  angle 

GAM7RI 

Y 

rad/sec 

Rate  of  change  of  elevation  flight 
path  angle 

HGC7F 

h 

ft 

Geode t i c  a  1 1 i tude 

PH  1  PD 

*P 

deg 

Angle  between  body  axis  and  platform 
ax!  s 
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Symbol 

Quant i ty 

Un  i  ts 

Remarks 

P177R 

O 

rad/ sec 

Inertial  angular  rate  of  body  x  axis 

PS  1  PD 

^P 

deg 

Angle  between  body  axis  and  platform 
axis 

Q177R 

q 

rad/ sec 

Inertial  angular  rate  of  body,  y  axis 

RG77N 

R 

9 

naut.  miles 

Approx,  range  from  launch  point 

R177R 

r 

rad/sec 

Inertial  angular  rate  of  body,  z  axis 

SIG7D 

a 

deg 

Horizontal  flight  path  angle 

SIG7R1 

a 

rad/sec 

Rate  of  change  of  horizontal  flight 
path  angle 

Tl'TPD 

0 

P 

deg 

Angle  between  body  axis  and  platform 
axis 

THTRD 

9 

r 

deg 

Pitch  angle  between  rotating  machinery 
axis  and  body  axis 

TIME 

t 

sec 

FI  ight  Time 

TIMES 

ts 

sec 

Stage  Time 

U777F 

u 

ft/sec 

Inertial  velocity  component  along  body 
x  axis 

V777F 

V 

ft/sec 

Inertial  velocity  component  along  body 
y  axis 

VA77F 

VA 

ft/sec 

A  i  r speed 

VG77F 

*  <J 

ft/sec 

Ground  referenced  speed 

W777F 

w 

f  t/sec 

Initial  velocity  component  along 
body  z  axis 

WTR7P 

w 

t 

1  bs 

Weight  of  body 

XG/7F 

X 

g 

ft 

Downrange 

XG77FI 

X 

g 

f  t/sec 

Rate  of  ch. nqp  of  downrange 

YG77F 

Y 

g 

ft 

C  ros  s  range 

YG77F! 

Y 

g 

t t/sec 

Rato  of  change  of  crossrange 

ZG77FI 

"g 

ft/sec 

Rate  of  change  of  altitude 
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Aerodynamic  Forces  and  Moments  -  SACS 


( 1 )  1  n  1 1 

ial  Print 

None 

(2)  Time 

History  Print 

Symbol  ( 

Quantity  Units 

Remarks 

CA 

CA 

Axial  force  coefficient 

CAVAH 

o 

> 

o 

1 

Ca  at  a  =  6  =  0° 

CM 

Cm 

Pitching  moment  coefficient 

CN 

CN 

Normal  force  coefficient 

CRM 

C| 

Rolling  moment  coefficient 

CY 

Cy 

Side  force  coefficient  (body  axis) 

CYM 

c 

n 

Yawing  moment  coefficient 

Engine  Thrust  and  Fuel  Flow  - 

TITS 

( 1 )  1  n  1 1  i 

ial  Print 

None 

( 2 )  Time 

History  Print 

Symbol  Quantity  Units 

Remarks 

Mt 

M-j-  ft- lbs 

Net  engine  pitch  moment 

NT 

Nj  ft-lbs 

Net  engine  yaw  moment 

TO) 

T;  lbs 

Actual  thrust  array 
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Vehicle 

Physical 

Character i st i 

ics  -  VPCS 

(1)  Initial  Print 

Symbol 

Quant  i  t^ 

f  Units 

Remarks 

AREFF 

S 

ft2 

Reference  area 

DIRFF 

d, 

ft 

Reference  length  for  longitudinal 
aerodynamics 

D2RFF 

d2 

ft 

Reference  length  for  lateral  and 
directional  aerodynamics 

XCGRF 

XC.G. 

REF  ft 

Longitudinal  position  of  center  of 
gravi ty 

(2)  Time  Histor> 

'  Print 

Symbol 

Quant i tv 

'  Units 

Remarks 

Af FOL-TR-  71-1 55 

Pan  .  '  ?  ; 

3.  AUTOPILOT  OUTPUT 

The  autopilot  printout  is  divided  wito  nine  caf.-‘>ories:  auxiliary 
computations,  maneuver  logic,  pitch  autopilot,  yaw  autopilot  roll 
autopilot,  throttle  autopilot,  brak"  autopilot,  control  response,  and 
status  indicators.  Depersdinc;  on  the  values  of  the  autopilot  output 
indicators,  (see  Suction  il  Para.  3r  pg.  2h)  any  combination  of  the  nine 
categories  can  appear  in  the  autopilot  output. 

a.  Auxiliary  Computat ions  Output  (AUXICP=!.! 

The  output  variables  here  are  X,  (ft),  YR  (ft),  ZR  (ft),  XR  (ft/  , 
YR  (ft/sec),  ZR  (f../sec),  ip  (deg/sec),  and  <£p  (deg/sec)  (see  Ref.  i, 
pg  226)  whose  column  des : g~at ions  are  XR,  VR,  ZR,  XRDl ,  YRD1  ,  ZRDI, 
PSIPDl,  and  PHIPD1  respectively. 

b.  Maneuver  Logic  Output  (MANLiJ>G=n 

The  output  va^’afales  h*rre  depend  on  whether  the  problem  phase  is 
glide  slope,  flare,  landing  roll  or  takeoff  roll. 

(1)  Glide  Slope  Output  (IAP=1) 

The  variables  01  .  .it  here  are  \Jr  (ft/se  :) ,  Vad(ft/sec)  ,  YR(deg), 

QR(lb/ft2),  CLR,  ad  (deg)  ,  LR(lbs),  URObs),  Td(lbs'l,  hGS(ft),  RR(ft), 

h  ( f  t  j  ,  h  (ft),  h  (ft),  h  (ft),  h  (ft),  and  <}>  (deg)  (see  J.  I, 
ea  e  ei  pa  pt  a 

Dgs  )60  &  161)  whose  column  designations  are  VE  ,  \IM,  OAMPR,  QR  f'LR, 
ALPDES,  LR,  OR,  TO,  HGS ,  RR ,  HEA ,  HE,  HET,  HPT,  and  PHIDES  respectively. 

(2)  Flare  Output  (IAP*2,3) 

The  variables  output  here  are  t^(sec),  th(sec),  xyp(ft),  h^fft)  , 

V  (ft/sec)  ,  Vh|,,;  ft/sec)  ,  XRf.(ft),  XRfr(ft/sec),  hRfr(ft),  hR.  (ft/sec), 

A  3 ( f t/sec  ,  A  (ft/sec?),  u.(deg),  T  (lbs),  and  c,(deg)  {,  C  (ft) 
xK  hK  c  d  dm 

(Ref  I,  g  1 73~  1 76)  *»nose  column  des  Ignat  ic  iS  art  TX,  TH,  XTD,  HTO,  VXTD, 
VHTO ,  Xrtf .  XRFO 1 ,  HRF ,  HRFDI,  AXR ,  AMR.  ALPDES,  TD,  PHIDES,  and  OH, 

'  *SPeCt  ively.  As  the  (  lore  an,  r->  the  hold  nxide,  only  ;  and  T 

a .  •  output. 
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(3)  Landing  Roll  Output  ( I  AP=4 ) 

The  variables  output  here  are  ISS,  ILR,  ICS.  IBS,  T|(sec),  and 
t^Csec)  (Reference  I,  pg  195)  whose  column  designations  are  ISS,  ILR, 
ICS,  IBS,  Tl,  and  TR,  respectively. 

(4)  Takeoff  Roll  (UP-5,  6) 

The>e  is  no  specific  output  for  the  ii.aueuver  logic  of  the  takeoff 
rol I  phase . 

c.  Pitch  Autopilot  Output  (PITCHP=1) 

The  variables  output  here  are  i^ldeg),  ag(deg)  ,  a  (&)  •  U&l  ■ 
ae^.(deg),  and  6^(deg)  (Reference  1,  pg  189)  whose  column  designations 
are  DELQN ,  ALPdAE,  ALPHdl ,  ALPDD1  ALPHET  ,  and  DELQDE .  respect i ve lv , 

d.  Yaw  Autopilot  lutout  (YAWAUP=1) 

The  variables  output  here  are  6  (deg),  3  (deg),  I  '  ,  8  ..  (deg), 

rn  c  ;^scc  j  6 1 

S^'deg),  ^(deg),  and  ^^(deg)  (Reference  i,  pg  193)  whose  column 
designs  ions  are  DELQR,  BETAD ,  BETAD1  ,  BET"  F.T ,  DElRDE,  PSIE,  and  PSIET, 
respect i ve ly . 

e.  Pol  Autopilot  Output  (R$LLAP«1 ) 

The  variables  outpu-  here  are  4'e(deg),  (deg)  ,  and  ^^^(deg) 

(Re  fere. -ce  1,  oc  195)  whose  column  designations  ar  PHIE,  PHItT,  and 
DELPDE,  respectively. 

f.  Throttle  Autopilot  Output  ( T  HRg>AP=  I  > 

The  '  ab  i  es  output  r  .ere  are  the  two  arrays  N^O)  and  T^(t> 
(Reference  1,  pgs  ’96-216)  whose  column  desig"  'ons  are  ND  (IN)  and 
TO  (IN),  f esnec t i ve !y . 
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g.  Brake  Autopilot  Output  (BRAKAP=l) 

The  basic  output  variable  here  is  the  braking  moment  array  Mg.(lb-ft). 
When  the  braking  signal  (|RS  =  1)  is  given  in  the  landing  roll,  arrays 

“m  (fir)  *  “tr!  (—2)  *  and  ^Ei  (fir)  (Reference  pg  221>  are 

also  output.  The  arr?v  output  names  are  MB ( 1 ) ,  OMEGATR(l),  QMEGATR9 1 ( I ) , 
and  OMEGATE(l),  respectively. 

h.  Control  Response  'C<t>NTRP=l) 

The  variables  output  here  are  5 •  ,  6*  <5.  '^£2.) 

q  '  sec  J  r  '  sec  '  p  sec ' 

6^(deg),  <S ^ ( deg )  ,  6^(deg)  and  the  arrays  N^.^  nnd  N^.^  (Reference  1, 
pg  224)  whose  column  designations  are  DELQD1,  DELRDi,  DELPD1,  DELQD, 

DELRD,  DELPD,  ND1,  ana  N,  respectively. 


i.  Status  Indicator  Output  ( I  NO  I CP= 1 ) 

The  problem  phase  indicator  IAP  and  the  indicator  arrays  I C  ( I )  and 
lu.  are  output  here  und-r  the  names  IAP,  I C  ( I )  and  IB(l),  respectively. 

D  I 

4.  LANDING  GEAR  OUTPUT 

Ihe  gears  must  be  staged  in  (see  pg  69)  for  gear  output  to  occur. 
The  kind  and  amount  of  gear  output  depends  on  the  indicator  INOLG. 

a.  INOLG  =  -I 

This  value  of  the  indicator  is  used  to  stage  in  the  gears  and  do 
all  calculations  required  to  obtain  strut  positions,  velocities,  and 
wheel  speeds  Gear  output  under  this  condition  is  as  foil  ,<s : 


V  triable 

Output  Name 

*yft) 

OELTA 

pK',b/ft2) 

P 

P^'Vt2) 

P2 

r,K(.bs) 

FT 
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Variable 

Output  Name 

Variable 

Output 

RK (  /seC  ^ 

SR 

FTRA(,bs> 

FTRA 

SFK(lbs) 

SF 

FTRB(,bs> 

FTRB 

Q^\  /sec1") 

AA 

FTRCObs) 

FTRC 

FCK2(lbs) 

FC2 

htx ( 1 b-ft) 

MTX 

UK 

MUVP 

mty (Ib-ft) 

MTY 

VGPTK^  /SCC^ 

VGPT 

MT2(lb-ft) 

MTZ 

FTRXK(,bs) 

FTRX 

FXM(lbs) 

FXM 

FTRYK(,b5) 

FTRY 

FyM(lbs) 

FYM 

FTR3tJ1bS* 

FTRZ 

F3M(ibs) 

FZM 

MAK(,b'ft) 

MA 

LM(lb-ft) 

LM 

V(,b’ft) 

MB 

MM(lb-ft) 

MM 

6k  (ft) 

DDELTA 

nm  (l b-f r ) 

NM 

5K(Ft/sec^) 

S02 

SK(Ft/sec) 

SD! 

sK(ft) 

S 

SK2 (Ft/sec2) 

S2D2 

SK2 (f l/sec) 

S20 1 

sK2(ft) 

$2 

u)TK(rad/sec2) 

OMETO) 

i.krK(rad/sec) 

OMET 
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All  variables  are  explicitly  defined  in  Reference  1  except  +  6^ 
and  +  6^  is  only  the  positive  value  of  6^  (Reference  1,  Eq.  199).  If 
6^  is  negative,  +  6^  is  given  a  zero  value;  +  6^  represents  the  actual 
tire  deflection.  is  the  sum  of  forces  resisting  main  strut  movement, 
that  is 


c  sc  -  P  4*F  -C§ 

SFK  AK  CK2  K2  K2 


K2 


-  c  s  + 

tK2L  k2 


-  F  (see  Eq  159  Ref.  !) 


b.  INDLG  =  2 

This  value  of  the  indicator  is  used  in  the  Stiff  Strut  Stage.  Since 

the  strut  positions,  velocities,  and  wheel  speeds  are  fixed  in  this  stage, 

only  the  following  output  is  made:  +(^K’  FTK’  ^K’  FTRXK'  FTRYK’  FTR3K’ 

V  SK’  FTRA’  FTRB’  FTRC’  Mrx’  MTZ’  FXM’  FYK’  F3M’  lm’  mm 

and  Nm- 


The  output  column  names  remain  unchanged  from  the  previous  case. 

Several  closing  comments  are  appropriate  as  regards  the  output.  As 
was  mentioned  on  pg  A8,  the  autopilot  calculations  are  only  updated  every 
time  interval  DElTS.  If  a  PRTMIN  is  reached,  an  autopilot  output  may  or 
may  not  occur  depending  on  how  PRTMIN  falls  within  the  DELTS  interval. 
Autopilot  output  will  always  appear  on  a  PRINT  condition.  However,  to 
make  the  autopilot  output  pertain  to  the  specific  data  in  the  SDF-2 
output,  the  PRINT  interval  must  be  the  same  as  DELTS  or  an  integer 
multiple  of  DELTS  . 
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SECTION  VI  I 
PROGRAM  USE 

The  amount  of  specific  data  needed  in  the  data  deck  depends  on  the 
problem  phase  to  be  studied.  TOLA  was  designed  for  conventional,  powered 
aircraft  and  therefore  has  certain  limitations  when  applied  to  uncon¬ 
ventional  vehicles.  These  limitations  are  also  closely  associated  with 
the  problem  phase  to  be  studied. 


GLIDE  SLOPE 


There  are  two  baste  glide  slope  assumptions:  The  glide  slope 
elevation  angle,  e^,  is  small;  the  vehicle  is  powered.  Both  are  valid 
assumptions  for  conventional,  powered  aircraft.  As  such,  the  powered 
vehicle  can  maintain  a  specific  glide  slope  position  arid  velocity  down 
the  glide  slope  subject  to  different  wind  conditions  and  aerodynamic 
changes  due  to  ground  effect.  An  unpowered  vehicle,  however,  does  not 
have  this  glide  slope  control  flexibility.  For  a  given  longitudinal 
trim,  the  glide  slope  angle  and  ground  speed  are  fixed  and  depend  on 
wind  conditions.  The  logic  for  glide  slope  control  of  such  a  vehicle  is 
n^t  in  TOLA.  At  best,  TOLA  can  be  used  to  confirm  unpowered  vehicle 
glide  slope  performance  for  specific  known  steady  state  conditions.  As 
such,  use  of  the  TOLA  glide  slope  phase  is  predominantly  limited  to 
conventional,  powered  vehicles. 


To  do  a  glide  slope  calculation  requires  specific  SDF-2  data  input 

and  specific  autopilot  data  input.  The  only  landing  gear  data  needed  is 

the  R  .  i nput . 
gR 

a.  Glide  Slope  SDF-2  Data  Input 

All  basic  SUl -2  data  input  (see  Section  II  para  I)  is  needed.  The 
aircraft  must  be  started  close  to  the  desired  glide  slope  position  and 
velocity  and  be  trimmed  approximately  for  the  desired  conditions.  This 
requires  inputs  of  the  proper  magnitude  for  the  following  SDF-2  variables: 
HGC7F,  THT8D,  GAM/D,  VG77F,  ATAB$I ,  and  RGR .  If  winds  are  involved,  the 
inputs  must  include  the  WTAb  tables. 


81 


AFFDL-TR-71-155 
Part  III 


b.  Glide  Slope  Autopilot  Data  Input 

All  engine  data  (see  Section  II  para  3a)  must  be  input.  The  throttle 
settings,  N,  should  be  close  to  that  required  for  trimmed  power.  No 
drag  chute  data  (see  Section  II  para  3b)  is  needed.  Only  ITO,  HF,  and  NF 
are  required  in  the  phase  begin  data  (see  Section  II  para  3c).  ITO  must 
have  a  zero  value  and  NF  s  juld  have  a  one  value  to  end  the  calculation 
at  the  HF  altitude.  The  other  phase  begin  input  data  is  not  needed 
unless  one  desires  the  calculation  to  go  into  the  flare  phase,  hold 
mode,  etc.  No  takeoff  data  is  needed  (see  Section  II  para  3d). 

All  glide  slope  maneuver  data  is  needed  (see  Section  I!  para  3e) . 

The  data  input  here  requiring  careful  scrutiny  is  ALPHDL,  DELEPS,  DELSIG, 
RFH,  PGS,  RFY,  and  PHIC.  ALPHDL  should  not  exceed  the  stall  conditions 
for  the  aircraft.  DELEPS  and  DELSIG  determine  the  accuracy  of  the  glide 
slope  position  control  and  should  be  approximately  0.1  degree.  To  see 
the  long  period  motion  after  a  perturbation  (i.e.,  wind  changes,  engine 
failures,  etc.),  PGS  and  PHIC  must  be  zero.  If  it  is  desired  to  control 
the  long  period  oscillations,  RFH,  PGS,  RFY,  and  PHIC  must  be  appro¬ 
priately  selected  in  sign  and  magnitude,  which  can  be  selected  easily 
with  a  knowledge  of  basic  aircraft  performance  and  the  aircraft  input 
data.  The  flare,  hold,  and  landing  roll  maneuver  data  (see  Section  II 
para  3f,  3g.  3h)  is  not  needed.  Engine  failure  in  the  glide  slope  is 
staged  through  the  variables  HI  and  H2  (see  Section  II  para  3i)  and 
their  associated  arrays  IHI  and  IH2.  No  bi ake  condition  data  (see 
Section  II  para  3j)  is  needed  for  the  glide  slope.  The  four  autopilots 
pitch,  yaw,  roll,  and  throttle  must  be  built  by  judicious  selections 
of  the  constants  associated  with  each. 

(1)  Pitch  Autopilot.  The  data  input  here  (see  Section  II 
para  3k)  is,  for  the  most  part,  self  explanatory.  That  data  input 
requiring  careful  examination  is  ALPDL ,  DELALA,  RFALPH,  and  PSH.  ALPDL 
is  used  to  detect  a  discontinuous  change  in  (for  example,  on  wind 
changes  and  transition  from  glide  slope  to  flare  logic)  thru  the  term, 
and  prevent  the  large  rate  from  entering  the  autopilot  (see  Appendix  I). 
The  ALPDL  value  must  be  larger  than  the  expected  in  the  flare. 
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Since  the  discontinuities  ot^  produce  very  large  rates,  it  is  safe  to 
pick  a  sizeable  value  for  ALPDL  (i.e.,  like  2  or  3  deg/sec)  and  be 
assured  that  legitimate  signals  are  not  limited.  DELALA  determines 
the  accuracy  of  the  angle  of  attack  control.  A  DELALA  value  of  0.01-0.02 
degree  is  sufficient  for  good  control.  The  constants  RFALPH  and  PSH 
determine  the  pitch  control  of  the  aircraft.  It  is  left  up  to  the  user  to 
choose  appropriate  values  for  the  particular  aircraft  system  being 
simulated.  Suffice  it  to  say,  that  values  can  be  easily  obtained  (with 
a  basic  knowledge  of  aircraft  pitch  performance  and  the  aircraft  input 
data)  that  control  the  aircraft  to  the  command. 

(2)  Yaw  Autopilot.  The  data  input  here  (see  Section  II  para  3 l) 
is  self  explanatory.  As  with  the  pitch  autopilot,  the  appropriate  values 
of  the  constants  RFB,  PSR,  RFPSI,  and  PSPS  I,  depend  on  the  aircraft  system 
and  are  to  be  determined  by  the  user. 

(3)  Roll  Autopilot.  Comments  on  the  data  input  here  (see 
Section  II  para  3m)  is  much  the  same  as  with  the  yaw  autopilot. 

(4)  Throttle  Autopilot.  The  data  input  here  (see  Section  II 
para  3n)  is  self  explanatory. 

No  brake  autopilot  data  input  (see  Section  II  para  3o)  is  needed 
for  the  glide  slope.  The  control  response  data  input  (Section  II  para  3p) 
is  neeed  for  the  glide  slope;  the  data  is  self  explanatory.  The 
initialization  input  data  (see  Section  II  para  3q)  and  the  autopilot 
output  indicators  (Section  II  para  3r)  are  self  explanatory.  No  specific 
staging  (see  Section  II  para  4)  Is  needed  in  the  glide  slope  except  the 
cards  at  the  end  of  the  data  deck  (see  Section  II  para  4g) . 

2.  FLARE 

The  flare  guidance  laws  develop  a  constant  acceleration  maneuver 
appropriately  limited  by  touchdown  conditions.  To  achieve  the  requested 
acceleration  vector  in  general  requires  two  control  degrees  of  freedom 
(i.e.,  angle  of  attack  and  power).  Upper  and  lower  bounds  exist  on  both 
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angle  of  attack  and  power,  so  a  situation  can  exist  where  the  requested 
acceleration  vector  is  outside  the  control  capability  cf  the  aircraft. 

In  such  a  case,  the  control  variable  magnitudes  used  are  those  that  give 
the  least  vector  error  in  the  requested  acceleration  vector.  The  additional 
input  (i.e.,  other  than  that  required  for  glide  slope)  for  the  flare 
is  as  follows: 

All  gear  input  data  (see  Section  II  para  2)  is  needed.  In  the  phase 
begin  input  (see  Section  II  para  3c)  NF  must  be  given  a  zero  value.  As 
previously  discussed  (Section  II  para  3c),  DELTAH  should  always  have  a 
nonzero,  positive  value  (i.e.,  like  1.  ft).  However,  if  a  decrab  or  tail 
down  constraint  maneuver  is  expected,  DELTAH  should  be  large  enough 
(note  DELTAH  is  approximately  the  distance  between  the  bottom  tire 
surface  of  the  main  gears  and  the  runway)  to  give  the  autopilot  sufficient 
time  to  perform  the  required  maneuver.  NLRI  should  have  the  value  one 
to  stop  the  calculation  on  impact. 

All  flare  maneuver  input  data  (see  Section  II  para  30  is  needed. 

That  data  requiring  careful  input  is  DA,  TL ,  and  TU.  For  an  accurate 
flare  search,  DA  should  be  0.01-0.02  degree,  etc.  DA  should  at  least 
bf  as  small  as  OELALA  in  the  pitch  autopilot.  TL  and  TU  determine  the 
thrust  range  allowed  during  the  flare  maneuver.  If  the  vehicle  is  not 
powered,  TL  and  TU  must  be  given  zero  values.  In  the  hold  maneuver  data 
(see  Section  II  para  3g) i  the  kill  engine  indicator  array,  KE ,  and  the 
tail  down  constraint  PM  are  needed.  Engine  failure  (see  Section  II 
para  3i)  in  the  flare  is  staged  on  the  runway  altitudes  HR1,  HR2,  and 
their  associated  arrays  IHR1  and  IHR2.  In  the  pitch  autopilot  input  (see 
Section  II  para  3k),  PSH2,  RFALP2,  and  DELQC2  must  be  added.  Note  that 
in  the  yaw  autopilot  data  input,  the  overcontrol  constants  PSR  and  PSPSI 
must  have  opposite  signs.  Two  stages  are  required  to  get  the  flare  to 
terminate  on  runway  impact:  stage  gears  into  program;  smooth  impact  stage 
(see  Section  II  paras  Aa  and  Ab) . 
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To  start  the  problem  in  the  flare  requires  the  data  input  of  HGC7F 
to  be  (Section  II  para  1)  less  than  the  flare  altitude  HF  (see  Section  II 
para  3c).  To  start  the  problem  in  the  hold  mode  requires  the  appropriate 
input  of  XRF  and  HRF  (Section  II  para  3c  and  Reference  1,  pg  182)  and  the 
input  of  ALPDES  and  TTD  (Section  II  para  3g) • 

3.  LANDING  ROLL 

The  vehicle  impact  section  of  the  landing  roll  problem  is  general 
in  that  it  can  simulate  the  ground  impact  of  a  rigid  main  structure 
with  up  to  five  independent  oleo  struts.  The  control  logic  for  the 
landing  roll,  however,  is  specialized  to  aerodynamic  lifting  vehicles 
which  land  horizontally.  Impact  of  unconventional  aircraft  and  vehicle 3 
which  do  not  iand  horizontally  can  be  determined,  however,  under  the 
restriction  that  all  control  variables  remain  fixed. 

a.  Horizontal  Landing 

If  a  drag  chute  is  employed,  the  chute  data  (Section  II  para  3b) 
must  be  input.  If  the  problem  starts  after  impact,  KP  must  have  a  one 
value  (Section  II  para  3c)  and  the  time  of  impact,  Tl,  must  be  given  an 
initial  value.  The  indicator  NLRI  must  have  a  zem  value.  The  runway 
stopping  conditions  VS,  XS,  and  TS  must  also  be  input. 


All  landing  roll  maneuver  data  (Section  II  para  3h)  must  be  input. 
Engine  failure  in  the  landing  roll  is  sequenced  through  the  variables 
TR1  and  TR2  (Section  II  para  3>)  and  their  associated  indicator  arrays 
ITR1  and  ITR2.  The  Drake  condition  input  data  (Section  II  para  3j )  is 
needed.  In  the  pitch  autopilot  data  (Section  I!  para  3k),  TST,  DELQF, 
and  DELFD1  must  be  input.  The  brake  autopilot  data  (Section  II  para  3o) 
is  needed.  Note  that  the  MBL  amy  should  always  contain  zeros  and  that 
DELTAW  must  be  small  enough  to  give  accurate  control  around  PD  (the 
relationship  between  DELTAW  and  PD  depends  on  the  shape  of  the  coefficient 
of  friction  -  percent  skid  table,  FTAB03  around  the  value  PD).  The  first 
three  stages  (Section  II  paras  La,  Lb,  and  Lc)  are  a  must.  The  spoiler 
aero  stage  is  needed  if  spoilers  are  used  in  the  landing  roll.  The  rest 
of  the  stages  can  be  included  at  the  users  descretion. 
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b.  Fixed  Control  Variables 

The  usefulness  of  this  option  for  unconventional  vehicles  (such  as 
STOL  and  VSTOL)  depends  upon  the  reality  of  the  fixed  control  variable 
assumption.  A  constant  body  oriented  thrust  vector  can  be  obtained  by 
input  of  Dl.FXP,  DLFYP,  and  DLFZP  in  the  SDF2  data  (see  pg  17)  along  with 
an  INDTFF  value  at  zero  (see  pg  52;  this  will  zero  the  thrust  received 
from  the  thrust  routine).  Fixed  pitch,  yaw,  and  roll  trim  can  be  obtained 
by  eliminating  all  input  of  aerodynamic  coefficients  associated  with 
elevator,  rudder,  and  aileron  defections  (note  this  automatically  zeros 
these  coefficient)  and  by  appropriate  input  of  ATAB10,  ATAB24,  ATAB38, 
ATAB51,  ATAB65,  and  ATAB80  (see  pgs  26-30).  The  aerodynamics  can  be 
completely  deleted  by  an  INDAER  indicator  value  of  zero  (s±e  pg  51) • 


4.  TAKEOFF  ROLL 


The  takeoff  roll  should  start  from  a  near  equilibrium  condition  for 
the  aircraft  strut  system  with  the  engines  at  the  takeoff  throttle 
setting.  The  takeoff  termination  altitude,  HS  (Section  II  para  3c  -  this 
is  normally  the  altitude  to  clear  a  50-ft  obstacle)  is  needed.  The  takeoff 
condition  data  (Section  II  para  3d)  is  needed.  Engine  failure  during 
takeoff  can  be  staged  on  the  variables  XRFI  and  XP.F2  (Section  II  para  3i) 
through  their  associated  arrays  IT1  and  IT2.  DELQTO  must  be  added  to  the 
pitch  autopilot  (Section  I  I  para  3k).  The  first  three  stages  (Section  II 
paras  4a,  4b  and  4c)  are  a  murt. 
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SECTION  VIII 
DECK  SETUP 

1.  DECK  STRUCTURE 

Running  the  TOLA  computer  program  requires  a  particular  deck  setup. 

The  deck  structure  is  presented  as  a  guide  only  in  determining  this  setup. 


An  end  of  record  (EOR)  card  is  a  7.  8.  9  punched  in  column  1  and  an  end 
of  job  (EOJ)  card  is  a  6,  7.  8,  9  punched  in  column  1. 
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2.  CONTROL  CAROS 

A.  Execute  TOLA  where  TOLA  is  a  binary  file  on  permanent  file 
TOLACP ,  Cyle  1. 


(Job  card) 

ATTACH, TOLA, TOLACP , I D-XXXXXXX ,CY«1 . 

TOLA. 

(end  of  record) 

[Data  cards] 

(end  of  job) 

B.  Ixecute  TOLA  where  TOLA  is  a  binary  file  on  permanent  file 
TOLACP,  Cycle  1;  generate  a  data  tape;  and  generate  a  CALCOMP  plot  tape 
from  the  data  tape.  The  plot  program  is  a  binary  file  on  permanent  file 
TOLAPLT,  Cycle  1. 

(Job  card) 

ATTACH, TOLA, TOLACP, I  D-XXXXXXX ,CY-1 . 

LABEL , TAPE  1 3  >  W,L  ,-PLTDATA .VSN-LXXXXX .  RING  IN 
TOLA. 

REWIND, TAPE13. 

REQUEST, TAPE7, HI ,  N , VSN-LXXXXX .  RING  IN 
ATTACH , PLTOLA .TOLAPLT , I D-XXXXXXX ,CY- 1 . 

PLTOLA. 

(end  of  record) 
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[Data  cards  for  TOLA] 

(end  of  record) 

[Data  cards  for  PLTOLA] 

(end  of  job) 

C.  If  the  TOLA  computer  program  is  on  MT,  replace  the  "ATTACH, TOLA, 
TOLACP, ID*XXXXXXX,CY“1 card  with  the  following  card  in  examples  a  and  b 
above : 


REQUEST, TOLA, MT,E,VSN=LXXXXX.  RING  OUT 
3.  CALCOMP  PLOTTING  INPUT 

The  following  data  is  required  in  order  to  generate  CALCOMP  plots 
by  the  PLOT  Tape  Generating  Program  (PLTOLA). 

A.  Data  nerated  on  file  TAPEI3  (disk  or  tape)  by  TOLA.  The  input 
required  by  TOLA  to  generate  data  on  file  TAPE  1 3  is  described  in  V.3.S., 
page  6A. 

B.  The  following  may  be  read  from  cards  on  the  input  file  using  the 
NAMELIST  feature  of  Fortran  Extended  with  the  group  name  of  INPUT.  If  the 
value  of  any  variable  is  the  same  as  its  nominal  value,  it  is  not 
necessary  to  read  it  as  input. 
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VARIABLE 

NOMINAL 

NAME 

VALUES 

VALUES 

NCASES 

1 

ISDFR 

1 

1 

0 

ISDF 

0 

0 

1 

ISTPR(i) 

1 

1 

0 

ISTPL(i) 

0 

0 

1 

IL 

0 

0 

l 

TFIRST 

0. 

TLAST 

0. 

PLTINT 

1 

K 

FCTR 

1.0 

XL 

7-2 

YL 

5.0 

DESCRIPTION 

Number  of  sets  of  data  or  cases  to 
be  plotted 

Rigid  body  data  is  stored  on  tape 
Rigid  body  data  is  not  stored  on 
tape 

Do  not  plot  rigid  body  data 
Plot  rigid  body  data 

Landing  gear  data  for  gear  i  is 
stored  on  tape 

Landing  gear  data  for  gear  i  is 
not  stored  on  tape 

Do  not  plot  landing  gear  data 
for  landing  gear  i . 

Do  plot  landing  gear  data  for 
landing  gear  i . 

Do  plot  lower  chamber  pressure 
and  2nd  piston  plots. 

Do  not  plot  lower  chamber  pressure 
and  2nd  piston  plots. 

Trajectory  time  to  begin  plotting 

Trajectory  time  to  stop  plotting 
If  both  TFIRST-TLAST-O. ,  the  entire 
time  history  on  tape  will  be 
plotted . 

Plot  every  Kth  point 

The  current  factor  all  coordinates 
are  multiplied  by.  That  is,  the 
plot  is  made  larger  or  smaller  if 
FCTR  is  greater  than  1.  or  less 
than  1.  For  example,  if  it  is 
desired  that  the  plots  to  be  25 % 
of  the  original  size,  let  FCTR*. 25 

Length  of  X-axis  of  plot  in  inches 

Length  of  Y-axis  of  plot  in  inches 
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T .  Some  examples  of  data  input  are  as  follows: 

(1)  Example  No.  1.  Plot  rigid  body  variables  and  landing  gear 
variables  for  gears  1,  3,  and  5.  Plot  every  point,  and  plot  entire  time 
history.  Assume  rigid  body  and  landing  gear  variables  for  gears  1,3, 
and  5  are  stored  on  TAPE13-  The  input  will  be  as  follows: 

$  INPUT  I SDF=1 , ISTPRH ,0 , 1 ,0, 1 , I STPL-1  ,0,1 ,0,1$ 

(2)  Example  No.  2.  Plot  landing  gear  variables  for  gear  No.  3* 
Plot  every  other  point  from  time  =  A.  to  time  =  10.  seconds.  Assume 
rigid  body  and  landing  gear  variables  for  gears  1,2, 3, A,  and  5  are  stored 
on  TAPE13  for  time  =  0  to  20.  seconds.  The  input  will  be  as  follows: 

$  I NPUT  I  STPL  (3)=*=  1  ,  TFIRST=A.,  TLAST=10.,  PLTINT=2$ 

(3)  Example  No.  3.  Plot  rigid  body  variables  and  landing  gear 
variables  for  gear  No.  5-  Plot  every  point  and  plot  entire  time  history. 
Assume  rigid  body  and  landing  gear  variables  for  gears  1,3,  and  5  are 
stored  on  tape  13-  Desire  that  the  size  of  graphs  to  be  50%  of  the 
original  size  where  the  original  size  of  the  X-axis  is  8  inches  and  the 
Y-axis  is  6  inches.  The  input  will  be  as  follows: 


$  I NPUT  I SDF=1 , 1 STPR=1 ,0,1 ,0, 1 , I STPL (5)-l ,XL=8. ,YL=6. ,FCTR«. 5$ 
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